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Revised  report  for  July  2004-June  2005 


Introduction 

Purpose  of  Research:  To  develop  theoretical  models  of  prion  protein  oligomeric 
aggregation,  strains,  clearance,  and  propagation  in  assistance  to  the  goal  of  the  overall 
Army  prion  project  of  developing  early  detection  schemes.  The  philosophy  is  that  a 
theoretical  understanding  of  these  different  elements  of  prion  science  can  provide 
potential  strategies  for  early  detection. 

Subject  of  Research  (referenced  to  original  approved  SOW): 

Theoretical  modeling  of  oligomeric  aggregation  of  prion  proteins  (Task  1  from  approved 
scope  of  work)  -  Pertinent  accomplishments  in  Apps.  B,C,D,E 

Theoretical  modeling  of  Prion  Strain  Phenomenology  (Task  II  from  approved  scope  of 
work)  -  Pertinent  accomplishments  in  App.  B 

Models  for  Prion  Clearance  (Task  III  of  approved  SOW)  -  Pertinent  accomplishments  in 
App.  E 

Models  for  Prion  Propagation  (Task  IV  of  approved  SOW)  -  Pertinent  accomplishments 
in  Apps.  B,D,E 

All  papers  with  full  citations  are  contained  in  the  appendix  list  at  the  end  of  this  revised 
report. 

Overview 

In  this  second  year  of  our  award,  we  have  made  significant  progress  on  several  areas  of 
theoretical  research  on  the  prion  problem,  particularly  on  the  structure  of  the  minimal 
prion  oligomer  and  the  potential  role  of  domain  swapping  in  that  structure,  the  role  of 
metal  binding  in  potentially  protecting  against  conformational  change,  and  the  role  of  the 
membrane  in  determining  the  kinetics  of  the  disease  time  course.  In  terms  of  our 
productivity,  we  have:  (i)  submitted  five  papers  for  publication,  one  of  which  have  been 
published,  one  accepted,  and  three  are  under  review;  (ii)  presented  seven  invited  lectures 
at  universities  or  meetings  and  three  contributed  talks  at  conferences;  (iii)  co-organized  a 
major  international  workshop  on  amyloid  diseases  (http://i2cam.org/i2camyloid);  (iv) 
received  national  press  coverage  for  our  research  on  prion  diseases  (see  the  attached  item 
in  Appendix  );  (v)  a  student  (Ms.  Jianping  Pan)  has  received  two  different  fellowships  to 
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support  her  research  in  this  area. 


In  the  coming  year  and  final  year  of  the  award,  building  upon  the  successes  emerging 
from  the  just  completed  sabbatical  leave  of  the  PI,  we  plan  to  explore  (i)  detailed 
molecular  models  for  strains  in  prion  disease  to  be  followed  by  kinetic  modeling  of 
strains;  (ii)  the  potential  role  of  the  proposed  domain  swapped  prion  trimer  model  in 
explaining  how  at  least  three  genetic  mutations  lead  to  disease;  (iii)  binding  of  other 
metal  ions  to  the  prion  protein  (notably  iron  and  sodium)  and  binding  of  copper  to  a 
different  site  away  from  the  octarepeat  regions  of  the  protein;  (iv)  exploration  of  a 
potential  role  of  the  octarepeats  in  converting  the  cellular  prion  protein  to  scrapies  form, 
in  formation  of  non-scrapies  based  oligomers,  and  in  providing  proteinase  resistance 
upon  copper  binding;  (v)  modeling  the  potential  role  of  the  neuronal  membrane  in 
mediating  fission  or  autocatalysis  of  the  scrapies  form  of  the  prion  protein. 


Key  Research  Accomplishments 

1.  Proposing  a  new  structure  for  the  minimal  infectious  unit  of  prion  disease:  the 
Domain  Swapped  Prion  Trimer  (DSTP).  [General  relevance  to  overall  goal  of 
development  of  an  oligomeric  aggregation  model  (Task  I  of  approved  SOW);  specific 
relevance  to  Task  2.d  -  structural  modeling  of  prion  strains.]  This  work  has  been 
accepted  for  publication  in  the  Hypothesis  section  of  FASEB  Journal  and  is  presented 
in  App.  B.  One  of  the  more  exciting  developments  in  the  study  of  prions  in  recent 
years  has  been  the  discovery  of  areal  aggregates  of  infectious  oligomers  in  electron 
microscopy  work  by  the  UCSF  group  (see  refs.  3,  4  of  App.  B).  In  Ref.  4  of  App.  B, 
the  UCSF  group  boldly  proposes  that  the  infectious  oligomers  are  trimers  of  left 
handed  (3-helices  with  a  similar  shape  and  size  to  the  bacterial  carbonic  anhydrase 
trimer  (ref.  5  of  App.  B),  and  that  these  oligomers  can  readily  stack  into  fibrils.  The 
difficulty  with  this  proposal  lies  in  understanding  what  stabilizes  the  proposed 
(3-helical  Prion  Trimer  (BPT):  while  the  carbonic  anhydrase  possesses  three  zinc  ions 
which  covalently  bond  to  histidines  on  different  monomers,  there  is  nothing 
analogous  for  the  BPT.  We  note  in  this  paper  that  domain  swapping,  implicated 
elsewhere  in  amyloid  formation  (see  Refs.  11-13  of  App.  B)  can  stabilize  the  trimer 
by  affording  hydrogen  bonding  in  the  central  region  (Fig.  1  of  App.  B).  Using  all 
atom  MD  simulations,  we  showed  that  this  structure  is  held  together  by  a  significant 
hydrogen  bonding  network  (Fig.  2  of  App.  B)  similar  to  the  serine  acyltransferase 
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trimer  (ref.  9  of  App.  B,  also  studied  as  a  control  in  our  paper),  compatible  with  the 
external  a-helices  of  the  BPT  model  (Fig.  3  of  App.  B),  is  remarkably  stable  to 
temperature  (which  blows  the  BPT  structure  apart — see  Fig.  4  of  App.  B),  and  may 
provide  insights  into  the  conformational  details  of  prion  strains  (p.  19  and  Fig.  5  of 
App.  B)  and  three  mutations  implicated  in  heritable  prion  disease  (discussion  on  p.  19 
of  App.  B).  We  discuss  these  latter  points  further  in  our  discussion  of  research  plans 
for  the  coming  year.  We  regard  the  DSTP  model  as  an  important  modification  to  the 
proposed  BPT  structure  for  the  infectious  prion  protein.  This  work  was  carried  out  in 
collaboration  with  S.  Yang,  H.  Levine,  and  J.  Onuchic  from  the  Center  for 
Theoretical  Biological  Physics  at  UC  San  Diego,  a  fruitful  collaboration  for  the  PI 
which  will  be  carried  forward  in  the  coming  year. 

2.  Comprehensive  computational  study  of  metal  binding  to  the  prion  protein  and 
demonstration  of  a  potential  role  of  copper  binding  in  inhibiting  conversion  of  the 
prion  protein.  [General  relevance  to  Task  l  of  approved  scope  of  work  -  metal 
binding  can  provide  potential  detection  mechanisms  via  ESR  and  SPECT,  and  can 
disrupt  aggregation  as  per  the  arguments  in  the  paper.]  [Note:  a  shortened  version 
of  this  work  discussing  only  the  influence  of  copper  on  inhibiting  aggregation 
appeared  in  Biophys.  J.  ]  In  this  work,  with  manuscript  submitted  to  Journal  of 
Molecular  Biology  and  presented  in  App.  C,  we  have  performed  a  fully  Quantum 
Mechanical  electronic  structure  study  of  the  binding  of  several  divalent  transition 
metal  ions  to  the  prion  protein  and  explored  the  potential  implications  of  such  metal¬ 
binding  to  the  normal  function  of  the  prion-protein  and  to  their  role  in  prion-diseases. 
We  find  that  metal  ions  bind  to  at  least  two  different  regions  of  the  prion  protein:  (i) 
the  highly  conserved  HGGGW  sequence  in  the  octarepeat  region  near  the  N-terminal 
and  (ii)  the  GGGTH  sequence  in  the  region  of  the  protein  that  is  believed  to  convert 
to  the  beta- sheet  form  in  course  of  the  prion  diseases  (For  location  of  the  sequences  in 
the  human  prion  protein  See  Fig.  1  of  the  manuscript  in  App.  C).  For  the  HGGGW 
sequence,  the  binding  energy  trend  found  in  our  studies  for  different  transition  metal 
ions,  reported  in  Table  2  of  the  manuscript  in  App.  C,  is  consistent  with  in  vitro 
measurements,  and  can  be  plausibly  explained  by  Jahn-Teller  and  electrostatic 
mechanisms  compatible  with  the  well  known  Irving-Williams  series  for  binding 
affinities  of  transition  metal  ions  with  polydentate  ligands.  The  coordination 
geometry  of  copper  ion  in  the  octapeptide  region  (See  Table  1  and  Figs.  2  and  3  of 
manuscript  in  App.  C),  agrees  well  with  Electron  Paramagnetic  Resonance  (EPR)  and 
X-ray  crystallography  measurements.  We  also  find  that  GGGTH  sequence  binds 
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much  more  strongly,  consistent  with  separate  experiments.  Two  of  our  findings  may 
have  particularly  important  implications  for  the  (normal)  cellular  function  and 
disease-causing  role  of  the  prion  protein:  (i)  We  find  that  metal-binding  affinity  and 
associated  peptide  conformations  are  strongly  dependent  on  protonation  states  of  the 
peptides,  implying  that  pH  changes  can  modulate  metal-binding  and  protein 
conformation  (See  Figs.  2  and  3  of  the  manuscript  in  App.  C).  (ii)  We  find  that  the 
metal-bound  conformation  is  incompatible  with  the  misfolded  disease  causing 
(infectious)  structure  proposed  by  Govaerts  et  al  (See  Ref.  16  and  Fig.  5  of  the 
manuscript  in  App.  C  for  comparison  details).  These  findings  suggest  that  copper 
needs  to  be  expelled  from  the  protein  before  the  beta-sheet  rich  misfolding  can  occur. 
This  is,  indeed,  consistent  with  experimental  findings  of  copper  depletion  in  infected 
and  diseased  tissues  (ref.  20  of  App.  C)  and  that  copper  suppresses  in  vitro  amyloid 
formation.  Overall,  our  studies  are  consistent  with  the  idea  that  metal  ions,  such  as 
copper,  play  a  complex  role  in  prion  diseases.  They  may  play  a  protective  role  in  the 
disease  by  creating  a  disease-incompatible  structure  and  thus  enhancing  the  barrier  to 
misfolding.  However,  the  structure  may  be  protease  resistant,  and  this  may  enhance 
disease  propagation.  This  work  has  been  submitted  to  the  Journal  of  Molecular 
Biology. 

3.  Detailed  modeling  of  in  vitro  yeast  prion  aggregation  [Of  general  relevance  to  Task 
I  of  original  SOW,  of  specific  relevance  to  Task  3.b  (developing  spatially  explicit 
kinetic  models),  Task  4.b  (development  of  kinetic  models  for  propagation),  and  as  a 
warmup  to  task  2.b,  developing  kinetic  models  of  prion  strains.]  The  importance  of 
yeast  prions,  in  the  field  of  prion  diseases,  arises  from  the  fact  that  several  key 
concepts  of  mammalian  prion  propagation  have  been  validated,  by  studying  relatively 
harmless  yeast  proteins  in  the  laboratory.  These  include  fission  and  exponential 
growth,  in  vitro  generation  of  infectious  agent,  and  existence  and  propagation  of 
strains.  Detailed  modeling  of  yeast-prion  kinetics  allows  one  to  highlight  the 
similarities  and  differences  between  yeast  and  mammalian  prions,  which  may  be 
important  for  controlling  the  disease.  In  this  work,  with  manuscript  submitted  to 
Physical  Review  E  and  presented  in  App.  D,  we  present  a  model  of  yeast  prion 
aggregation  that  can  account  for  almost  all  aspects  of  in  vitro  kinetic  data  that  show 
(i)  Linear  aggregates  which  grow  by  monomer  addition,  (ii)  A  modest  dependence  of 
the  lag-time  on  monomer  concentration  varying  nearly  as  an  inverse  square-root 
function  (See  Fig.  2  in  the  manuscript  presented  in  App.  D),  (iii)  Sigmoidal  growth  of 
aggregates  that  is  nearly  exponential,  and  (iv)  A  distribution  of  aggregate  sizes  with  a 
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mode  that  is  large  (10s- 100s  of  monomers).  We  have  contrasted  our  simulation  work 
done  via  kinetic  equations  (Figs.  3, 4, 5, 6  in  the  manuscript  presented  in  App.  D)  and 
by  stochastic  algorithms  (  Figs.  7,8,9  in  the  manuscript  in  App.  D)  with  alternative 
models  of  aggregation  that  do  not  include  fission  and  hence  fail  to  obtain  the 
exponential  growth.  We  have  also  shown  that  fission  leads  to  a  distribution  of 
aggregate  sizes  that  is  Weibullian  (See  Fig.  4  and  Fig.  9  in  the  manuscript  presented 
in  App.  D),  which  may  be  tested  in  future  experiments.  In  comparing  the  yeast-prion 
kinetics  to  the  mammalian  case,  we  find  two  important  differences.  The  first  is  that 
mammalian  prions  may  form  two-dimensional  areal  aggregates  rather  than  one¬ 
dimensional  ones.  A  more  significant  difference  seems  to  be  the  role  of  membranes  in 
fissioning  and  exponential  growth,  which  is  not  necessary  in  the  yeast  case,  allowing 
spontaneous  fission  to  take  place  in  vitro.  It  is  clearly  important  to  further  investigate 
the  role  of  membrane  in  fissioning  of  mammalian  prions. 

4.  Demonstration  of  the  necessity  of  the  membrane  in  exponential  in  vivo  growth  of 
prion  infectivity.  [General  relevance  to  Task  I  of  SOW  (infectious  oligomer 
aggregation  requires  membrane  binding  to  spread  disease),  general  relevance  to 
Task  III  (clearance  requires  membrane  binding),  general  relevance  to  Task  IV 
(disease  propagation  with  attendant  exponential  growth  requires  membrane 
binding)]  [Note:  this  work  was  published  in  Biophys.  J.  ]  This  work  has  been 
submitted  as  a  technical  comment  to  Science  and  is  contained  in  App.  E.  We  respond 
here  to  a  beautiful  paper  by  Chesebro  et  al.  (ref.  1  of  App.  E)  showing  that  membrane 
anchoring  is  necessary  for  cell  toxicity/clinical  symptoms  in  mammalian  prion 
disease  by  working  with  transgenic  mice  lacking  the  sequence  necessary  for  attaching 
the  GPI  anchor  of  the  prion  protein.  At  the  same  time,  the  animals  remain  infectious! 
We  point  out  several  issues  in  our  comment:  (i)  Data  from  Chesebro  et  al.  showing 
the  time  course  of  prion  infectivity  in  the  transgenic  mice  is  very  well  fit  by  a 
quadratic-in-time  dependence  over  the  full  600  days  of  observation  (!)  (Fig.  1  and  Eq. 
1  of  App.  E),  consistent  with  simple  elongation  kinetics  described  by  Ferrone  (refs.  3, 
4  of  App.  E)  without  (a)  monomer  depletion  or  (b)  fission/exponential  growth.  The 
significance  of  this  result  is  that  the  membrane  anchoring  must  play  a  crucial  role  in 
generating  the  exponential  growth  characteristic  of  prion  disease,  (ii)  We  extend  the 
argument  to  infer  that  the  fundamental  elongation  rate  for  wild  type  aggregates  (with 
membrane  anchors)  is  likely  to  be  significantly  slower  than  in  the  anchorless  case  and 
attribute  this  to  what  is  most  probably  a  slower  clearance  rate  for  the  anchorless  PrPc 
proteins  (Eq.  2  and  discussion  on  p.  76  of  App.  E).  This  result  helps  to  explain  why 
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mice  heterozygous  for  anchor-full  and  anchorless  PrPc  display  faster  incubation 
times  at  a  given  dose  than  wild  type  mice. 


5.  Review  of  amyloidogenic  proteins  from  a  materials  perspective.  In  this  paper, 
presented  in  App.  F  and  published  in  the  June  2005  issue  of  the  Materials  Research 
Bulletin,  we  have  put  forward  a  review  of  the  proteins  implicated  in  amyloid  diseases 
from  a  materials  science  perspective.  While  this  does  not  represent  new  research,  it 
does  represent  a  significant  intellectual  effort  aimed  at  communicating  our  research 
interests  to  a  much  broader  scientific  audience. 


Research  Directions  for  Coming  Year 

1.  Correlation  of  relative  mutant  hydrophobicity  with  incubation  time  and  DSTP 
stability.  We  have  initiated  discussions  with  Glenn  Telling  and  Karah  Nazor  of 
the  University  of  Kentucky  about  understanding  their  Tg  mice  data  for  mutations 
at  P 101  using  our  DSTP  model.  They  find  in  mice  overexpressing  the  PrP  protein 
that  mutation  at  P101  (implicated,  e.g.,  in  the  P102L  mutation  leading  to  GSS  in 
humans)  to  residues  more  hydrophobic  than  P  (notably  L  and  I)  induces 
spontaneous  prion  disease  with  incubation  times  roughly  correlating  with 
hydrophobicity.  Within  our  DSTP  model,  this  may  be  potentially  be  rationalized 
by  the  following  observations:  (i)  for  PrPc,  substituting  a  relatively  more 
hydrophobic  residue  for  PI 02  will  partially  destabilize  the  structure  since  this  is  in 
a  solvent  exposed  region;  (ii)  for  our  DSTP  PrPSc  model,  the  substitution  can  have 
two  effects,  (a)  generically  rendering  the  proline  containing  hinge  region  of  our 
DSTP  more  flexible  and  thus  potentially  lowering  the  barrier  to  conversion  by 
increasing  the  search  rate  of  the  protein,  (b)  affording  solvent  exclusion  in  the 
core  region  and  thus  stabilizing  the  DSTP  relative  to  the  WT  form.  In  the  latter 
case,  a  simple  two  state  model  (which  is  admittedly  an  oversimplification) 
displays  a  linear  relationship  between  the  change  in  the  free  energy  difference 
between  cellular  and  scrapies  form  of  the  protein  and  the  lowering  of  the  barrier. 
Our  goal  here  is  to  carry  out  all  atom  MD  simulations  on  both  cellular  and  DSTP 
forms  of  the  prion  protein  and  compare  to  the  incubation  time  data  of  Telling  and 
Nazor  to  assess  the  molecular  basis  for  the  correlation  with  relative 
hydrophobicity.  This  will  be  carried  out  in  collaboration  with  S.  Yang,  H. 
Levine,  and  J.  Onuchic  at  UCSD. 
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2.  Modeling  of  strain  conformations  and  kinetics  within  the  domain  swapping 
picture.  As  noted  in  App.  B,  the  introduction  of  domain  swapping  opens  up 
intriguing  conformational  possibilities  for  the  explanation  of  strains  in  prion 
disease,  the  simplest  of  which  is  that  the  number  of  monomers  in  a  domain 
swapped  oligomer  encodes  strain.  If  this  is  so,  strains  will  be  able  to  breed  true 
by  either  (i)  oriented  areal  aggregation  on  a  membrane  surface  or  (ii)  by 
interneuronal  templating  across  (most  likely)  the  synaptic  cleft.  In  addition  to  this 
simple  idea,  there  are  three  other  ideas  worth  exploring  within  this  general 
picture:  (a)  Strain  conformation  is  encoded  within  the  'splay’  of  the  left  handed  13- 
helices — the  helices  of  carbonic  anhydrase  are  for  example  rather  straight,  while 
those  for  serine  acyltransferase  are  splayed  inward  from  C-terminus  to  N- 
terminus.  This  idea  could  account  for  why  some  strains  have  associated  amyloid 
plaques  and  others  don’t,  since  splay  will  mitigate  against  stacking  due  to  the 
epitaxial  mismatching  between  the  helical  ends,  (b)  Strain  conformation  is 
encoded  in  the  cross  sectional  shape  of  the  left  handed  (3-helix.  While  most  of  the 
20  non-redundant  left  handed  (3-helices  found  on  the  pdb  possess  triangular  cross 
section,  one  has  two  sides  (the  collagen  binding  domain  of  yersinia  adhesion 
YadA,  pdb  code  1P9H).  It  is  worth  exploring  whether  such  a  two  sided  helix 
might  be  possible  for  the  prion  protein,  (c)  Domain  swapping  occurs  with  the 
major  loop  of  the  prion.  Cohen  and  Prusiner  (Ref.  19  of  App.  B)  have  earlier 
suggested  that  different  domain  swaps  might  correspond  to  different  strains.  In 
this  light  it  is  interesting  to  consider  the  possibility  of  domain  swapping  via  the 
major  loop  of  the  BPT  structure.  This  would  be  ineffective  for  strains  which  can 
take  the  'miniprion’  PrPSc  106  form,  but  is  worth  exploring.  We  will  first  explore 
possibility  (a)  and  examine  in  particular  the  stability  of  different  dimers  and 
tetramers.  Assuming  we  find  stability  for  the  different  possible  strain 
conformations,  we  will  develop  stochastic  kinetic  models  for  both  areal 
aggregation  and  interneuronal  templating  to  assess  the  feasibility  of  this  simple 
picture  describing  time  course  kinetics. 

3.  Possible  docking  of  octarepeats  with  left  handed  /?- helices  and  implications  for 
cellular-to-scrapies  conversion.  Zahn  has  shown  that  the  octarepeat  region  of 
the  prion  protein  with  copper  bound  may  represent  a  new  structural  motif  (Ref.  2 1 
of  App.  C).  Intriguingly,  three  octarepeats  yield  a  structure  resembling  two  sides 
of  an  equilateral  triangle  with  cross  sectional  area  approximately  that  of  the 
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standard  left  handed  (3-helix.  This  suggests  a  potential  role  for  docking  of  the 
octarepeat  on  the  left  handed  (3-helix  of  a  scrapies  form  prion  which  could  be 
important  in  templating  the  cellular  form  of  the  prion  for  conversion.  We  intend 
to  study  this  docking  possibility  with  all  atom  molecular  dynamics  and 
subsequently  explore  the  implications  for  the  conversion  process.  We  note  that 
while  octarepeats  are  not  necessary  for  infection,  they  do  enhance  incubation  time 
relative  to  transgenic  animals  lacking  them,  and  additional  (>4)  octarepeats 
induce  heritable  disease.  Hence,  the  presence  of  octarepeats  clearly  enhances  the 
tendency  towards  conversion  to  scrapies  form. 

4.  Metal  ion  binding  and  structural  modulation  of  the  prion  protein.  Besides  the 
Cu,  Ni,  Mn,  Zn  ions  already  studied,  we  will  look  at  binding  of:  (i)  Fe  to  the  priori 
protein.  This  is  motivated  by  work  by  Singh  and  collaborators1  who  have  shown 
that  the  prion  protein  co-transports  across  epithelial  cells  with  ferritin,  (ii)  Na  to 
the  prion  protein.  This  is  motivated  by  the  work  of  Nishina  et  al.  (Ref.  48  of 
App.  C)  who  have  shown  that  in  solution  both  Na  and  Cu  can  confer  enhanced 
proteinase  K  resistance  to  the  truncated  PrP  27-30  scrapies  form.  Since  nominally 
Na  is  monovalent  and  larger  than  Cu,  the  binding  mechanism  if  any  is  mysterious 
and  we  will  examine  potential  binding  sites  and  their  implication  for  PK 
resistance.  In  addition,  we  will  explore  how  Cu  binds  to  the  Hill  residue  as 
claimed  in  Ref.  12  of  App.  C. 

5.  Modeling  of  membrane  mediated  fission  and  oligomeric  autocatalysis.  Our 
work  in  App.  E  strongly  suggests  that  the  membrane  plays  a  role  in  engendering 
exponential  growth  of  prion  infectivity.  This  can  either  arise  from  fission  of 
aggregates,  or  from  autocatalysis.  The  former  case  clearly  applies  for  in  vitro 
yeast  prion.  We  propose  to  explore  two  potential  mechanisms  for  membrane 
mediated  fission,  both  deriving  from  Ginsburg-Landau  modeling  of  the  membrane 
bound  prion  aggregates.  In  the  first,  we  will  examine  whether  curvature  near  the 
edges  of  the  synaptic  cleft  can  break  aggregates.  These  areal  aggregates  are 
presumably  less  strongly  bound  than  individual  oligomers.  A  separate  possibility 
is  that  membrane  undulations  can  break  up  rafts.  Our  colleague  Prof.  Atul  Parikh 
at  UC  Davis,  studying  the  raft  like  structures  which  can  bind  other  GPI  proteins 
has  come  to  this  conclusion.  We  can  initiate  some  combined  theory/experiment 
studies  with  Prof.  Parikh  on  this  matter.  In  addition,  assuming  there  is  some 
population  of  "soft  oligomers”  from  the  cellular  form  of  the  prion  protein,  there 
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can  be  autocatalytic  interneuronal  templating  of  scrapies  form,  and  we  will 
attempt  to  simulate  this. 


Additional  Research  Results 

Last  year  we  discussed  a  study  of  the  membrane  insertion  properties  of  the  prion  protein, 
particularly  because  the  known  transmembrane  segments  of  the  prion  proteins  contain  a 
mutation  (A  117V)  implicated  in  GSS.  We  carried  this  out  with  the  segment  of  prion 
protein  from  residue  90-153,  which  contains  both  known  transmembrane  domain 
(residues  1 1 1-134),  the  primary  (3-helical  region  of  the  BPT  and  DSTP  models  discussed 
above,  and  of  course,  residue  117.  We  did  study  this  using  the  membrane  insertion 
program  used  earlier  by  our  group  for  a  study  of  the  insertion  of  A|3  peptides.  However, 
we  found  nothing  using  the  program  which  was  surprising,  unlike  the  A(3.  In  particular, 
(i)  whether  initiated  in  the  membrane  or  outside  the  membrane,  the  transmembrane 
sequence  reached  the  same  a-helical  structure,  which  this  particular  code  is  biased 
towards  producing;  (ii)  the  A117V  mutation  has  no  significant  effect  within  the 
simulation. 

We  initiated  a  study  of  the  possibility  of  using  Cu  or  Mn  isotopes  for  either  PET  or 
SPECT  studies  for  early  prion  detection,  but  put  that  aside  when  the  DSTP  study  opened 
up  multiple  research  directions. 

Reportable  Outcomes:  Awards,  Activities,  Presentations  and  Publications 

1.  Awards,  Activities,  and  Presentations  (Appendix  A) 

2.  Structure  of  Infectious  Prions:  Stabilization  by  Domain  Swapping  (Appendix 
B)  Published  as  S.  Yang,  H.  Levine,  J.  Onuchic,  D.l.  Cox,  FASEB  J  ]9,  1778- 
1782(2005) 

3.  An  ab  initio  Electronic  Structure  Study  of  Transition  Metal  Binding  to  the 
Cellular  Prion  Protein  (Appendix  C)  Shorter  version  published  as  D.L.  Cox,  J. 
Pan,  R.R.P.  Singh,  "A  mechanism  for  copper  inhibition  of  infectious  prion 
conversion,”  Biophys.  J.  2L.L11-L13  (2006) 

4.  One  Dimensional  Model  of  Yeast  Prion  Aggregation  (Appendix  D  )  Published 
as  K.C.  Kunes,  D.L.  Cox,  R.R.P.  Singh,  Phys.  Rev.  E72,  051915  (2005) 

5.  Comment  on  “Anchorless  Prion  Protein  Results  in  Infectious  Amyloid 
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Disease  Without  Clinical  Scrapie”  (Appendix  E)  Published  as  D.L.  Cox, 
R.R.P.  Singh,  S.  Yang,  "Prion  Disease:  Exponential  Growth  requires 
Membrane  Binding”,  Biophys.  J.  90,  L77-79  (2006) 

6.  The  Materials  Science  of  Protein  Aggregation  (Appendix  F)  Published  as  D.L. 
Cox,  H.  Lashuel,  K.Y.  Lee,  and  R.R.P.  Singh,  Materials  Res.  Bull.  30,  79-84 
(2005) 

7.  Yahoo!  News  Search  Results  for  Cox  Singh  "mad  cow"  (Appendix  G) 


Conclusions 

In  summary,  this  year  we  have  (i)  produced  a  new  model  for  the  minimal  prion  oligomer; 
(ii)  modeled  the  affinity  trend  and  structure  of  copper  binding  to  the  prion  protein  and 
shown  that  the  conversion  to  left  handed  p-helical  structure  of  the  scrapies  form  requires 
expulsion  of  the  nonoctarepeat  copper  consistent  with  post  mortem  studies,  cultured 
neuron  studies,  and  in  vitro  amyloid  growth  studies;  (iii)  in  vitro  yeast  prion  kinetics  can 
be  well  described  by  a  simple  stochastic  growth  model  with  fission;  (iv)  the  membrane  is 
critical  for  exponential  growth  of  prion  infectivity  and  without  membrane  anchoring 
prion  aggregate  elongation  is  more  rapid.  We  propose  a  number  of  studies  following 
from  this  work  towards  (i)  understanding  the  conformational  nature  and  kinetics  of 
strains;  (ii)  understanding  how  mutations  to  PI 02, 105  can  induce  heritable  disease;  (iii) 
how  unusual  structure  of  the  octarepeat  region  can  play  a  role  in  templating  conversion  of 
cellular  to  scrapies  form  of  the  prion  disease;  (iv)  how  Fe  and  Na  ions  might  bind  to  the 
prion  protein,  and  how  copper  may  bind  to  HI  1 1;  (v)  how  the  membrane  can  play  a  role 
in  engendering  exponential  growth  of  infectivity. 
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APPENDIX  A:  ACTIVITIES,  HONORS,  AND  PRESENTATIONS  ASSOCIATED 
WITH  THIS  RESEARCH  AWARD 

Honors 

J.  Pan  I"CAM  Junior  Fellow  (supported  collaborative  visit  with  Prof.  E.  Artacho, 
Cambridge  University,  Jan-Feb.  2005) 

J.  Pan  UC  Davis  Summer  Graduate  Research  Fellow  (supports  one  month  of  summer 
research,  Summer  2005) 

Activities 

D.L.  Cox,  R.R.P.  Singh,  co-organizers  of  First  International  ICAM  Workshop  on 

PROTEIN  AGGREGATION  AND  AMYLOID  FORMATION  IN 
SYSTEMIC  AND  NEURODEGENERA TIVE  DISEASE: 

PHYSICAL,  MOLECULAR  AND  BIOLOGICAL  APPROACHES,”  Lausanne, 
Switzerland,  July  16-19  2005.  (http://i2cam.org/i2camyloid) 

Invited  Presentations 

D.L.  Cox: 

•  "Molecular  level  modeling  of  the  prion  protein,  ”  Brooks  group  seminar,  Scripps 
Research  Institute,  June  2005. 

•  "Consider  a  spherical  mad  cow:  a  physical  look  at  prions Condensed  matter 
theory  seminar,  Department  of  Physics,  The  Ohio  State  University,  April  2005. 

•  "Consider  a  spherical  mad  cow:  a  physical  look  at  prions,"  Condensed  Matter 
Journal  Club,  Department  of  Physics,  UC  San  Diego,  March  2005. 

•  "Physical  Modeling  of  Amyloid  Diseases,"  invited  lecture,  Symposium  on 
Frontiers  in  Computational  Biology,  Rice  University,  December  2004. 

•  "Physical  Modeling  of  Amyloid  Diseases,"  seminar,  Center  for  Theoretical 
Biological  Physics,  UC  San  Diego,  November  2004. 


R.R.P.  Singh: 

•  “ Prion  Diseases ”,  Seminar  at  Jawaharlal  Nehru  University,  Delhi,  India,  July 
2004. 

•  “The  mysterious  prion  protein ”,  Mathematics  Department  Seminar,  University  of 
California  at  Davis,  November  2004. 

J.  Pan: 

•  "Theoretical  study  of  metal  binding  to  the  cellular  prion  protein,’’  Seminar, 
Biophysics  Graduate  Group,  UC  Davis,  November  2004. 
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Contributed  Presentations: 


D.L.  Cox.  "Theoretical  investigation  of  metal  binding  affinities  in  the  prion  protein,” 
Keystone  Symposium  on  Transmissible  Spongiform  Encephalopathies,  Snowbird,  Utah, 
January  2005. 

K.  Kunes,  Kinetic  and  Stochastic  Models  of  ID  yeast  "prions" ,  American  Physical 
Society  March  Meeting,  Los  Angeles,  March  2005. 

(http  ://meeting  s .  aps  .org/Meeting/MAR05/E  vent/2 1 907) 

J.  Pan,  Theoretical  study  of  metal  binding  in  the  Prion  protein,  American  Physical 
Society  March  Meeting,  Los  Angeles,  March  2005. 

(http  ://meeting  s .  aps  .org/Meeting/MAR05/E  vent/2 1 906) 
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Structure  of  infectious  prions:  stabilization 
by  domain  swapping 

Sichun  Yang,*  Herbert  Levine,*  Jose  N.  Onuchic,*  and  Daniel  L.  Cox*’^’1 

*Center  for  Theoretical  Biological  Physics,  University  of  California  San  Diego,  La  Jolla,  California, 
USA;  and  ^Department  of  Physics,  University  of  California,  Davis,  California,  USA 


ABSTRACT  A  candidate  structure  for  the  minimal 
prion  infectious  unit  is  a  recendy  discovered  protein 
oligomer  modeled  as  a  (3-helical  prion  trimer  (BPT); 
BPTs  can  stack  to  form  cross-p  fibrils  and  may  provide 
insight  into  protein  aggregates  of  other  amyloid  dis¬ 
eases.  However,  the  BPT  lacks  a  clear  intermonomer 
binding  mechanism.  Here  we  propose  an  alternative 
domain-swapped  trimeric  prion  (DSTP)  model  and 
show  with  molecular  dynamics  (MD)  that  the  DSTP  has 
more  favorable  intermonomer  hydrogen  bonding  and 
proline  dihedral  strain  energy  than  the  BPT.  This  new 
structural  proposal  may  be  tested  by  lysine  and  N 
terminus  fluorescent  resonance  energy  transfer  (FRET) 
either  direcdy  on  recombinant  prion  protein  amyloid 
aggregates  or  on  synthetic  constructs  that  contain  the 
proline/lysine-rich  hinge  region  critical  for  domains  to 
swap.  In  addition,  the  domain  swapping  may  provide  1 ) 
intrinsic  entanglement,  which  can  contribute  to  the 
remarkable  temperature  stability  of  the  infectious 
prion  structure  and  help  explain  the  absence  of  PrPSc 
monomers,  2)  insight  into  why  specific  prolines  are 
potentially  relevant  to  three  inherited  forms  of  prion 
disease,  and  3)  a  simple  explanation  of  prion  strains 
assuming  the  strain  is  encoded  in  the  monomer  number 
of  the  oligomers. — Yang,  S.,  Levine,  H.,  Onuchic,  J.  N., 
Cox,  D.  L.  Structure  of  infectious  prions:  stabilization 
by  domain  swapping.  FASEB  J.  19,  1778-1782  (2005) 

Key  Words:  prion  •  domain  swapping  ■  (3  -helix  ■  amyloid 


MODEL  (5-HELICAL  PRION  TRIMER 

The  protein-only  model  for  prion  diseases  remains 
the  leading  candidate  to  explain  the  infectious  form  of 
the  disease  (1),  and  this  has  been  significantly  strength¬ 
ened  by  the  successful  production  of  disease  from 
synthetic  prion  protein  aggregates  (2).  Clearly,  the 
understanding  of  prion  disease  pathology,  propagation 
within  the  body,  and  therapeutic  approaches  will  be 
greatly  enhanced  by  a  detailed  knowledge  of  the  struc¬ 
ture  of  the  minimal  infectious  protein  unit. 

The  discovery  of  areal  aggregates  of  infectious  prion 
protein  oligomers  in  purified  brain  extracts  from  mice 
has  inspired  a  leading  candidate  model  for  the  minimal 
infectious  unit,  a  trimer  in  which  large  portions  of  the 
N  terminus  are  converted  to  left-handed  p-helices  (3, 


4).  We  call  this  the  |3-helical  prion  trimer  (BPT)  model. 
Figure  1 A  shows  a  top  view  of  the  BPT  in  the  (3-helical 
region.  This  BPT  displays  a  3-fold  symmetry  consistent 
with  maps  of  the  electron  density  difference  between  the 
PrP  27-30  oligomers  and  the  “miniprion”  PrPSc106  oli¬ 
gomers.  The  BPT  can  easily  be  stacked  to  fonn  filaments, 
which  can  then  be  wound  into  amyloid  fibrils.  The  model 
trimer  also  shares  detailed  structural  features  with  the 
bacterial  protein  carbonic  anhydrase  trimer  (5).  Inspired 
in  part  by  this  model,  others  have  explored  the  possibility 
that  the  [3-helix  plays  a  role  in  other  amyloid  structures 
such  as  aggregates  of  polyglutamine  peptides  (6),  (3-amy¬ 
loid  peptides  (7),  and  yeast  prion-like  proteins  (8). 

However,  it  is  not  clear  in  the  BPT  model  what  holds 
the  trimer  together.  The  carbonic  anhydrase  trimer, 
which  was  the  template  for  the  BPT  model  is  in  fact 
held  together  by  intermonomeric  covalent  bonding  of 
zinc  ions  (5),  which  are  not  present  in  PrPSc.  It  was 
suggested  that  a  hydrogen  bonding  network  is  respon¬ 
sible  for  this  binding  (4) ,  as  in  the  case  of  the  bacterial 
(3-helical  trimer  serine  acyltransferase  (SAT)  (9),  but 
no  effort  was  made  to  test  this  idea  quantitatively. 
Moreover,  the  yellow  loops  of  the  BPT  model  (which 
are  highly  conserved  (10))  shown  in  Fig.  1  A,  B  and 
expanded  in  Fig.  1 C  contain  two  prolines,  and  prolines 
tend  to  disfavor  this  type  of  bending. 


HYPOTHESIS:  A  DOMAIN-SWAPPED  PRION 
TRIMER  IS  THE  MINIMAL  INFECTIOUS  UNIT 

We  hypothesize  that  a  stable  oligomeric  structure  for 
the  infectious  prion  is  the  domain-swapped  trimeric 
prion  (DSTP)  of  Fig.  Iff,  which  can  increase  hydrogen 
bonding  and  reduce  elastic  energy  relative  to  the  BPT 
structure  proposed  earlier.  Further,  we  conjecture  that 
the  DSTP  structure  may  be  relevant  to  1)  the  unusual 
temperature  stability  of  the  prion  protein  and  the  lack 
of  observation  of  PrPSc  monomers  (both  due  to  the 
entanglement  of  the  DSTP  structure),  2)  the  kinetics  of 
at  least  three  inherited  forms  of  human  prion  disease, 
and  3)  the  encoding  of  prion  strains  in  conformation. 


1  Correspondence:  CTBP,  UC  San  Diego,  MC  0374,  9500 
Gilman  Dr.,  La  Jolla,  CA  92093-0374,  USA.  E-mail:  cox 
@  physics  .ucdavis.edu 
doi:  10.1096/fj.05-4067hyp 
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Figure  1.  A)  Proposed  fi-helical  structure  of  the  truncated 
PrPSc106  prion  trimer  (BPT)  taken  from  ref  4  (residues 
G89-F140).  B)  Structural  model  of  domain-swapped  trimeric 
prion  (DSTP)  of  the  present  paper.  C)  The  proposed  hinge 
region  expanded  from  the  yellow  loops  in  panels  A.  B. 

DOMAIN  SWAPPING:  MOTIVATION  AND 
RESULTS 

Concern  with  stability  of  the  BPT  model  led  us  to 
modify  it  by  domain  swapping  (11-13),  especially  since 
prolines  are  present  in  the  hinge  or  bridging  regions  of 
several  proteins  that  domain  swap  (14).  For  example,  in 
pl3sucl,  the  strain  on  two  prolines  in  the  hinge  mod¬ 
ulates  the  equilibrium  between  monomers  and  domain- 
swapped  dimers  (15,  16). 

We  constructed  our  DSTP  model  of  Fig.  IB  from  the 
proposed  BPT  structure  (we  thank  C.  Govaerts  for 
providing  his  coordinate  hie)  by  reorienting  the  mono¬ 
mers  to  place  the  loop  regions  near  the  trimer  center 
and  swapping  residues  89-99  to  a  neighboring  (3-helix 
in  a  cyclical  fashion,  with  the  BPT  loops  serving  as 
hinges.  This  DSTP  model  straightens  these  loops,  adds 
more  possible  hydrogen  bonds  in  the  loop  region,  and 
preserves  the  putative  monomer-monomer  spacing  and 
3-fold  symmetry  of  the  electron  density  maps  (3).  We 
ran  all-atom,  explicit  solvent  MD  (see  Appendix  for 
details)  on  energy  minimized  DSTP  and  BPT  structures 
using  AMBER8  (http://amber.scripps.edu).  As  a  con¬ 
trol,  we  applied  the  same  protocol  to  the  serine  acyl- 
transferase  (SAT)  (3-helical  trimer,  which  has  no 
strained  prolines  and  for  which  there  is  direct  crystal¬ 
lographic  evidence  for  intermonomer  hydrogen  bond- 
mediated  cohesion  (9).  Wee  also  carried  out  two  other 
simulation  tests:  1)  we  examined  the  stability  of  the 
DSTP  model  with  the  inclusion  of  the  C  terminus 
a-helices  proposed  in  ref  4,  and  2)  we  carried  out  three 
distinct  1-ns  runs  for  both  the  BPT  and  DSTP  structures 
to  make  sure  there  is  no  special  circumstance  with 
regard  to  individual  simulations. 


We  have  concluded  that  the  DSTP  is  significantly  more 
stable  than  the  BPT  and  therefore  a  more  promising 
candidate  structure.  First,  we  have  enumerated  direct  and 
water-bridged  hydrogen  bonds  between  monomers  in  die 
BPT  model,  our  DSTP  model,  and  the  SAT  structure  (see 
Appendix),  and  carried  this  out  for  three  separate  MD 
runs  on  each  structure.  We  show  the  direct  hydrogen 
bond  counts  for  three  different  1-ns  runs  of  the  BPT  and 
DSTP  structures  in  Fig.  2.  For  the  BPT,  we  find  2  ±  1, 
5  ±  1,  and  2  ±  1  direct  hydrogen  bonds  for  the  three  runs 
(averages  over  the  hist  100-ps  of  simulation  time),  and  at 
the  last  simulation  time  0,  0,  and  1  water-mediated  hydro¬ 
gen  bonds.  For  the  DSTP,  we  find  for  all  three  runs  9  ±  1 
direct  hydrogen  bonds  for  the  hist  100-ps  of  simulation 
time  and  2,  1,  and  1  water-mediated  hydrogen  bonds  in 
the  last  simulation  snapshot.  In  the  one  BPT  case  with  a 
relatively  high  number  of  direct  hydrogen  bonds,  the  N 
terminus  from  one  monomer  has  contacted  with  the  N 
terminus  region  of  another.  In  general,  the  BPT  hydro¬ 
gen  bonds  form  around  Q89  only,  while  the  DSTP 
hydrogen  bonds  form  in  the  hinge  region.  Note  that  Fig. 
2  shows  direct  hydrogen  bonds  formed  in  the  monomer- 
monomer  interface  for  the  DSTP  to  genetically  increase 
with  time  while  those  from  the  BPT  generically  decrease 
with  MD  time.  In  contrast  to  the  BPT,  the  simulated 
(3-helical  portion  (F140-A237)  of  the  SAT  trimer  has  11 
direct  and  3  water-bridged  hydrogen  bonds  (after  two 
stage  equilibration  with  no  restraints) . 

Second,  we  remark  that  the  domain  swapping  also 
relaxes  the  P101  and  P104  dihedral  energy  by  ~1.0 
kcal/mol  after  1-ns  of  MD.  Hence,  the  domain  swap¬ 
ping  helps  relieve  elastic  strain  as  anticipated. 

Third,  we  have  verified  that  the  DSTP  structure  for 
the  PrPs<  1  06  can  accommodate  the  C  terminus  a-h  di¬ 
ces  without  affecting  stability.  This  is  a  nontrivial  point, 
since  a  rigid  rotation  of  the  monomers  within  the  BPT 
is  not  possible  due  to  steric  constraints.  As  shown  in  Fig. 
3,  we  find  that  by  “peeling  off’  one  triangular  edge 
(residues  133-140)  of  the  lower  (C  terminus)  end  of 
each  (3-helix  to  form  a  random  coil  segment  that  we  can 
link  to  the  a-helices  that  begin  at  residue  177,  with  no 
resultant  change  in  direct  hydrogen  bonding  structure. 

Finally,  high  temperature  MD  provides  direct  evi¬ 
dence  for  BPT  instability;  at  500K,  the  DSTP  remains 
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Figure  2.  MD  simulations  show  that  the  nonswapped  BPT  (A) 
has  less  direct  hydrogen  bonds  in  the  interface  region  than 
the  domain-swapped  DSTP  (B) .  In  both  cases  we  have  simu¬ 
lated  residues  G89-F140.  Different  colors  correspond  to  three 
different  1  ns  MD  runs  in  each  case  initialized  at  different 
initial  temperatures  (see  Appendix) . 
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Figure  3.  A)  Model  for  DSTP  PrPSc106  (with  residues  G141- 
H176  deleted)  and  a-helices  included.  Residues  133-140  have 
been  removed  from  the  lower  rung  of  the  (3-helices  for  each 
monomer,  converted  to  random  coil,  and  attached  to  residue 
177  of  the  a-helices.  B)  Number  of  direct  hydrogen  bonds  for 
the  DSTP  structure  with  intact  a-helices.  The  a-helices  have  no 
effect  on  the  stabilization  induced  by  the  domain  swapping. 


compact  and  entangled,  while  the  BPT  monomer- 
monomer  separation  nearly  doubles.  Principal  compo¬ 
nent  analysis,  which  identifies  the  amplitude  and  char¬ 
acter  of  the  dominant  fluctuations,  shows  similar  BPT 
unbinding  tendencies  at  300K  as  in  the  500K  simula¬ 
tion,  but  no  unbinding  tendencies  for  the  DSTP.  Figure 
4  shows  the  radius  of  gyration  in  the  500K  runs  as  a 
function  of  simulation  time.  Within  the  1-ns  window  of 
simulation,  the  BPT  radius  appears  to  keep  growing, 
reaching  about  twice  that  of  the  apparently  more  stable 
DSTP  structure.  This  stabilization  through  entangle¬ 
ment  may  be  relevant  to  the  unusual  temperature 
stability  of  infectious  prion  proteins,  although  we  can¬ 
not  rule  out  that  such  stabilization  is  related  to  larger 
scale  aggregates  of  prion  oligomers.  The  entanglement 
may  also  be  relevant  to  the  inability  to  isolate  mono¬ 
meric  forms  of  PrPSc,  since  the  scrapies  form  is  stabi¬ 
lized  in  this  picture  only  in  aggregate. 


POSSIBLE  RELEVANCE  TO  INHERITED 
DISEASES  AND  PRION  STRAINS 

Finally,  we  discuss  the  possible  relevance  of  the  DSTP 
model  to  inherited  diseases  and  strains.  First,  proline 
point  mutations  yield  heritable  Creutzfeldt-Jakob  dis¬ 
ease  (P105T  for  humans,  P104T  for  mice)  or  Gerst- 
mann-Straussler-Scheinker  disease  (P102L,  P105L  for 


humans,  P101L,  P104L  for  mice;  see  chapter  14,  ref  1). 
Transgenic  mice  with  P101I  have  even  greater  disease 
susceptibility  whereas  those  with  PlOlYhave  somewhat 
less  susceptibility  than  PI  OIL  (K.  Nazor  and  G.  Telling, 
private  communication).  We  note  that  the  centrality  of 
these  sites  in  the  DSTP  hinges  affords  possible  burial  of 
the  relatively  hydrophobic  I,  L,  T,  and  Y  side  chains  and 
a  potential  speed  up  of  trimer  formation  through 
increased  hinge  flexibility.  Moreover,  the  insertion  of 
relatively  hydrophobic  residues  into  the  unstructured 
region  of  the  normal  PrPc  protein  will  raise  the  energy 
of  this  structure.  Indeed,  it  is  already  known  that  the 
P102L  mutation  is  marginally  less  stable  than  the 
wild-type  protein  (see  p.  694,  ref  1  ).  A  decrease  in  the 
stability  of  the  cellular  protein  form  relative  to  wild-type 
monomers  coupled  to  an  increase  in  stability  of  the 
scrapies  form  (associated  with  more  effective  hydro- 
phobic  burial  of  the  mutant  residue)  should  lead  to 
more  rapid  conversion  kinetics  via  ramping  of  the 
overall  free  energy  surface  slope  in  favor  of  the  scrapies 
form. 

In  Fig.  5  we  illustrate  two  ways  in  which  strains  might 
breed  true  if  the  monomer  content  within  an  oligomer 
encodes  the  strain  information.  With  this  assumption, 
strains  can  breed  true  either  by  1)  oligomer  conversion 
and  crystallization  on  membranes  (Fig.  5A),  with  local 
shape  matching  favoring  homogeneous  growth  in  a 
manner  analogous  to  the  oriented  aggregation  of  nano¬ 
particles  (17),  or  2)  via  interneuronal,  epitaxially  tem- 
plated  conversion  across  the  synaptic  cleft  or  interneu¬ 
ronal  gaps  (Fig.  5 B).  The  plausibility  of  the  latter 
conjecture  rests  on  the  length  of  unstructured  N  ter¬ 
minus  PrPc,  which  is  ~100  residues,  or  ~30  nm 
stretched  out,  easily  sufficient  to  span  the  interneuro¬ 
nal  gaps  of  10-20  nm.  Both  of  these  models  would 
enjoy  enhanced  conversion  kinetics  in  the  presence  of 
unconverted  or  “soft”  oligomers  of  PrPc  (18);  the 
interneuronal  mechanism  in  this  case  would  allow  for 
exponential  growth  via  oligomeric  autocatalytic  conver¬ 
sion.  We  note  that  a  different  domain  swapping  model 
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Figure  4.  The  entanglement  due  to  domain  swapping  intro¬ 
duces  additional  stability  for  the  DSTP  relative  to  the  BPT  as 
revealed  by  high  temperature  (T=500K)  simulations,  where 
the  radius  of  gyration  of  the  BPT  increases  more  rapidly  than 
that  of  the  DSTP,  which  reflects  unbinding  of  the  BPT. 
Principal  component  analysis  at  300K  reveals  a  similar  un¬ 
binding  tendency. 
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Figure  5.  Schematic  models  for  strains  breeding  true  under 
the  assumption  that  strain  is  encoded  in  monomer  number. 
For  simplicity,  we  assume  one  strain  is  encoded  as  trimers, 
one  as  tetramers.  A)  Areal  aggregation.  Given  local  size  or 
orientation  mismatches  between  trimers  (upper  image)  and 
tetramers  (lower  image),  energetics  will,  for  example,  disfa¬ 
vor  growth  of  trimers  off  of  tetramer  seeds  (lower  image).  B) 
Interneuronal  templating.  Assuming  that  seeds  template  con¬ 
version  across  interneuronal  gaps  or  the  synaptic  cleft,  trim¬ 
ers  will  seed  trimer  growth  (upper  image) ,  tetramers  will  seed 
tetramer  growth  (middle  image),  and  there  will  be  no  cross- 
seeding  between  trimers  and  tetramers. 

for  strains  has  been  conjectured  previously  (19);  sepa¬ 
rately,  domain-swapped  prion  dimers  have  been  crystal¬ 
lized  (20).  The  latter  result  is  consistent  with  prion 
oligomerization  via  intermolecular  disulfide  bonding  as 
has  been  shown  by  redox  induced  prion  hbrillization  in 
vitro  (21),  although  only  intramolecular  disulfide 
bonding  has  so  far  been  observed  in  vivo  (22).  Other 
specific  domain  swapping  mechanisms  than  those  we 
have  discussed  here  may  play  a  critical  role  in  amyloid 
fibril  formation  for  prions  as  has  been  observed  for 
other  proteins  (24-26).  For  the  PrP  27-30  scrapies 
form,  the  large  loop  (residues  E145-R163)  of  ref  4) 
made  from  some  of  the  deleted  residues  of  the 
PrPSc106  structure  could  serve  as  a  hinge  region  for 
alternative  domain-swapped  structures. 


PROPOSED  EXPERIMENTAL  TESTS 

While  the  DSTP  structure  may  ultimately  be  testable 
direcdy  on  recombinant  infectious  prion  oligomers,  as 


have  been  synthesized  by  Prusiner  and  collaborators 
recently  (2),  we  first  propose  testing  an  analog  system 
by  FRET  labeling.  Poly-Q  peptides  are  of  considerable 
interest  because  of  the  polyglutamine  diseases  (espe¬ 
cially  Huntington’s).  The  “critical”  length  near  36 
glutamine  repeats  for  disease  is  suggestive  of  a  left- 
handed  (3-helix  structure,  as  discussed  elsewhere  (4,  6, 
27).  The  key  experiment  we  propose  is  to  construct 
synthetic  peptides  of  the  form  Q12(KPSKPK)Qm, 
where  m  =  24-36.  The  goal  is  to  label  the  N  terminus 
Q’s  in  half  the  peptides  with  a  donor,  the  K’s  with  a 
different  donor,  and  in  the  other  half  put  in  N  termi¬ 
nus  Q  acceptors  and  K  acceptors.  If  the  argument  about 
stress  relaxation  and  hydrogen  bond  stabilization  works 
for  the  prion,  it  should  also  work  here.  There  should  be 
intense  K-K  FRET  given  intermonomer  separations  of 
the  side  chains  by  ~4  A  if  the  domain  swapping  arises. 
There  will  also  be  weak  N  terminus-Q-  N  terminus-Q 
FRET.  Presumably  there  could  be  an  N  terminus-Q-to- 
N  terminus-Q  FRET  signal  arising  with  time  if  the 
trimers  are  first  formed  and  then  they  fibrillize,  since 
the  corresponding  distances  would  get  reduced  by 
~60%  case  assuming  stacking  of  the  [3-helices.  If  this 
works  on  the  synthetic  peptide  it  should  be  attempted 
on  recombinant  prion  protein  prior  to  in  vitro  aggre¬ 
gation  as  per  ref  2.  We  would  anticipate  little  interfer¬ 
ing  contributions  from  other  lysines  in  the  prion  pro¬ 
tein  sequence  that  are  far  from  the  hinge  region. 


SUMMARY 

We  propose  a  domain-swapped  trimer  model  for  the 
minimal  infectious  prion  oligomer,  and  note  that  the 
domain  swapping  stabilizes  the  oligomer  relative  to  the 
earlier  proposed  model  by  reducing  stress  in  proline 
containing  loops,  increasing  intermonomer  hydrogen 
bonding,  and  promoting  entanglement.  We  conjecture 
that  domain  swapping  may  play  a  role  in  explaining  the 
temperature  stability  of  the  scrapies  form  of  the  prion 
protein,  the  absence  of  PrPSc  monomers,  inherited 
prion  disease  phenotypes  associated  with  mutating  pro¬ 
lines  in  hinges  of  the  DSTP,  and  that  conformation  may 
be  encoded  in  the  monomer  number  in  a  given  oli¬ 
gomer.  Finally,  we  propose  experimental  tests  based  on 
FRET  labeling  of  principally  lysine  residues  near  the 
crucial  pralines,  which  should  produce  intense  FRET 
response  if  the  domain  swapping  is  present.  pj| 


Appendix 

The  simulations  were  carried  out  with  the  AMBER8  molecu¬ 
lar  dynamics  package  (28)  with  the  AMBER  parm99  force 
field  and  explicit  TIP3P  waters,  following  these  procedures.  1 ) 
Before  we  start  a  molecular  dynamics  simulation,  we  perform 
an  energy  minimization  of  2000  steps  to  partially  relax  the 
entire  molecular  system.  2)  We  perform  a  two-stage  equilibra¬ 
tion  to  further  relax  the  protein  and  the  surrounding  solvent. 
In  the  first  stage,  we  start  the  system  from  a  low  temperature 
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of  approximately  100K  and  gradually  heat  up  to  300K  over 
20-ps  of  simulation  time.  We  perform  this  stage  of  equilibra¬ 
tion  with  the  volume  held  constant.  In  the  second  stage,  we 
equilibrate  the  system  using  pressure  and  temperature  con¬ 
trol  to  adjust  the  density  of  water  to  experimental  values.  We 
have  equilibrated  the  system  for  a  total  of  40-ps.  During  the 
first  two  steps,  the  backbone  heavy  atoms  of  the  (3-helical  (and 
also  the  a-helical  if  applicable)  portions  are  positionally 
restrained  using  a  harmonic  potential.  The  hinge  loops  and 
the  regions  linking  the  (3-helical  and  a-helical  domains  are 
allowed  to  freely  move.  To  achieve  the  three  different  MD 
runs  for  the  BPT  and  DSTP  structures,  we  start  with  different 
initial  temperatures  that  seed  a  different  run  through  a 
different  initial  thermal  distribution  of  atom  velocities.  Fi¬ 
nally,  we  carry  out  a  1-ns  production  run  at  constant  pressure 
and  temperature  with  no  positional  restraints.  For  the  high 
temperature  (T=500K)  simulations,  the  production  runs  are 
performed  at  constant  volume  with  a  reduced  time  step  of 
integration  to  prevent  any  possible  system  blowup. 

Hydrogen  bonds  are  calculated  based  on  the  following 
criteria:  1)  direct  hydrogen  bonds  in  interface  are  counted 
with  a  donor-acceptor  distance  of  3.5  A  and  a  hydrogen- 
donor-acceptor  angle  of  60°,  and  2)  water-bridged  hydrogen 
bonding  in  interface  is  counted  if  a  water  molecule  couples  by 
criterion  1 )  to  an  atom  from  each  of  two  monomers. 
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APPENDIX  C 


A  Mechanism  for  Copper  Inhibition  of  Infectious  Prion  Conversion 
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ABSTRACT  We  employ  ab  initio  electronic  structure  calculations  to  obtain  two  structural  models  for  copper  bound  in  the 
strongest  binding  site  of  the  noninfectious  form  of  the  prion  protein.  The  models  are  compatible  with  available  experimental 
constraints  from  electron  spin  resonance  data.  The  bending  of  the  peptide  backbone  attendant  with  the  copper  binding  is  not 
compatible  with  the  requisite  straight  /3-strand  backbone  structure  for  the  same  sequence  contained  in  two  recently  proposed 
models  of  the  prion  protein  structure  in  its  infectious  form.  We  hypothesize  that  copper  binding  at  this  site  is  protective  against 
conversion  to  the  infectious  form,  discuss  experimental  data  that  appear  to  support  and  conflict  with  our  hypothesis,  and 
propose  tests  using  recombinant  prion  protein,  genetically  modified  cultured  neurons,  and  transgenic  mice. 
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The  functional  role  of  the  normal  cellular  prion  protein 
(PrPc)  may  well  be  related  to  copper  binding,  potentially 
protecting  against  oxidative  damage  in  the  synaptic  region 
by  sequestering  divalent  copper  (1).  This  suggestion  is  sup¬ 
ported  by  data  for  transgenic  knockout  mice  devoid  of  the 
gene  for  expressing  PrPc;  these  mice  appear  to  suffer  late 
stage  oxidative  degradation  in  the  neuronal  synaptic  regions 
where  surface  bound  prion  protein  is  preferentially  concen¬ 
trated  in  wild-type  mice  (2,3).  The  strongest  copper  binding 
site  is  in  the  protein  region  that  converts  to  /3-sheet  structure 
in  the  infectious  (PrPSc)  form  (1).  We  use  electronic  structure 
calculations  to  study  two  possible  geometries  for  the  bound 
copper-PrPc  complex;  we  find  that  these  geometries  are  in¬ 
compatible  with  recently  proposed  models  (4,5)  for  PiPSl 
oligomers,  and  we  thus  hypothesize  a  mechanism  for  inhi¬ 
bition  of  PrPc-to-PrPSc  conversion  via  copper  binding. 

Copper  binds  principally  to  mammalian  PrPc  at  octarepeat 
sites  of  highly  conserved  form  (though  not  number)  in  a  now 
well-understood  pyramidal  geometry  with  binding  to  the 
peptide  backbone,  histidine  side  chains,  and  axial  waters  (6). 
These  octarepeats  (residues  60-91  in  humans)  are  not  es¬ 
sential  to  PrPSc,  which  remains  infectious  even  after  protein¬ 
ase  exposure  that  leaves  residues  92-230  intact.  One  strong 
copper  binding  site  is  present  in  the  PrPSc  core  region,  con¬ 
taining  the  sequence  92-96  GGGTH  for  humans.  Electron 
spin  resonance  data  suggest  there  is  binding  to  the  H96  (H95 
in  mice)  side  chain  and  the  G94  amide  group  (1).  It  is  not 
known  whether  the  primary  peptide  coordination  is  to  four 
nitrogens  (NNNN)  or  to  three  nitrogens  and  an  oxygen 
(NNON),  although  the  former  structure  has  been  conjectured 
to  be  more  likely  (1). 

We  have  studied  this  problem  computationally  with  the 
SIESTA  local  orbital-based  density  functional  theory  code 
using  conjugate  gradient-based  energy  minimization  to  ex¬ 
amine  possible  GGGTH-Cu(II)  geometries  (7,8).  We  have 
considered  the  NNNN  structure  and  an  NNON  structure 
compatible  with  the  ESR  data  (1)  and  our  results  are  shown 


in  the  upper  two  panels  of  Fig.  1  for  the  mouse  prion.  We 
built  initial  candidate  structures  of  the  form  Ac-Cu(H20)6 
(GGGTH)-NH2  using  ChemSite  Pro  (ChemSW,  Fairfield, 
CA)  and  VMD  (9).  For  each  geometry,  the  lowest  energy 
sampled  provides  our  model  bound  copper-prion  complex. 
We  used  double-zeta  (DZ)  basis  sets  for  light  atoms  (H,  C,  N, 
O)  and  double-zeta  polarized  orbitals  for  copper.  We  used 
the  Troullier-Martins  norm-conserving  pseudopotential.  We 
employed  the  Perdew-Becke-Emzerhof-based  generalized 
gradient  approximation  exchange  correlation  energy  func¬ 
tional.  Our  energy  cutoff  for  matrix  element  integration  was 
120  Rydberg.  All  calculations  were  carried  out  in  a  periodic 
unit  cell  size:  30.7  X  30.7  X  30.7  A.  Geometry  minimization 
was  carried  out  to  a  force  tolerance  of  32  pN. 

Our  main  result  is  inferred  from  comparing  the  two  pos¬ 
sible  GGGTH-Cu(II)  geometries  to  recently  proposed 
/3-helical  PrPSc  trimer  models  (4,5).  We  see  that  the  backbone 
bending  induced  by  copper  binding  is  not  compatible  with 
formation  of  /8-strands  in  the  left-handed  helices.  (It  is  un¬ 
clear  if  another  recently  proposed  model  has  the  same  char¬ 
acteristics  (10).)  Hence,  copper  in  the  nonoctarepeat  sequence 
GGGTH  can  protect  against  PrPc-to-PrPSe  conversion;  re¬ 
moval  of  this  copper  by  pH  reduction  associated  with  synaptic 
fluctuations  or  endocytosis  may  be  a  key  step  in  the  conversion 
pathway.  This  picture  is  supported  by  the  observation  that  only 
the  GGTH  sequence  is  required  for  binding,  and  eutherian 
prion  sequences  reveal  high  conservation  of  either  GGTH  or 
the  very  similar  GGSH  sequence  in  the  converting  region  (11). 

Removal  of  this  copper  by  pH  reduction  may  be  a  key  step 
in  the  PrPc-to-PrPSe  conversion  pathway  because  it  is  known 
that  acidic  conditions  favor  detachment  of  the  copper  (12). 
We  note  two  routes  for  such  a  pH  reduction  with  accom¬ 
panying  conversion  of  PrPc.  First,  upon  endocytosis  into  the 
cell,  the  pH  is  reduced.  Assuming  that  some  PrPSc  is  endo- 
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FIGURE  1  Potential  copper  binding  motifs  in  the  converting 
region  of  the  normal  (PrPc)  mouse  prion  protein,  which  are 
consistent  with  ESR  data  (1)  are  shown  in  panels  a  and  b.  The 
corresponding  copper-free  stretch  of  the  left-handed  /3-helix 
model  of  the  infectious  (PrPSc)  protein  from  Govaerts  et  al.  (4) 
is  shown  in  panel  c. 

cytosed  as  well  as  PrPc,  then  conversion  can  take  place  within 
the  endocytosed  compartments  after  the  copper  is  expelled. 
Second,  given  that  PrPc  is  preferentially  concentrated  near 
the  synapse,  then  from  signal  bearing  vesicles  released  near 
or  at  the  synapse  (13)  can  reduce  the  pH  locally  and  facilitate 
copper  expulsion. 

Three  experimental  results  support  our  proposed  mecha¬ 
nism  indirectly:  1),  Postmortem  studies  show  copper  deple¬ 
tion  in  regions  of  infection,  with  80-90%  reduction  of  copper 
content  in  the  corresponding  prion  proteins  (14).  This  result 
is  consistent  with  our  mechanism,  since  this  GGGTH  binds 
Cu(II)  first  with  potential  cooperative  enhancement  of  sub¬ 
sequent  binding  to  other  sites.  If  infection  promotes  removal 
of  copper  from  this  site,  then  it  is  more  likely  that  copper  will 
not  be  retained  at  the  other  sites  assuming  near-equilibrium 
conditions.  2),  Copper  uptake  in  infected  cultured  neurons  is 
suppressed  ~  10-fold  compared  to  control  cells  (15).  Again, 
this  is  consistent  with  the  hypothesis,  although  this  requires 
an  understanding  of  how  the  prion,  which  is  easily  shown 
to  contain  <1%  of  cellular  copper,  can  shut  down  copper 
transport.  3),  Copper  in  solution  inhibits  in  vitro  growth  of 
amyloid  fibrils  (16,17).  This  is  a  more  direct  corollary  to  our 
hypothesized  protective  role  of  copper. 

Other  experiments  provide  mixed  support  for  our  mech¬ 
anism.  The  copper  chelator  D-penicillamine  inoculated  in 
vivo  delays  disease  onset,  which  is  at  apparent  odds  with  our 
hypothesis  (18).  The  chelator  cuprizone,  however,  is  known 
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to  induce  a  spongiosis  (vacuolation  of  brain  tissue)  very 
similar  to  that  of  prion  disease  although  it  also  induces 
demyelination  of  neurons  and  no  transmissibility  (19,20). 
Moreover,  it  is  also  known  that  copper  confers  proteinase 
resistance  to  the  cellular  prion  protein  in  the  weaker  binding 
octarepeat  regions  (21).  This  is  significant  for  two  reasons. 
First,  proteinase  resistance  is  one  of  the  generic  hallmarks  of 
the  infectious  form  of  the  disease  (although  not  all  infectious 
prion  protein  is  in  fact  proteinase  resistant).  Second,  the 
octarepeat  deletion  does  not  remove  infectivity  from  PrPSc. 
However,  in  transgenic  mice  without  octarepeats,  disease 
incubation  upon  direct  intercerebral  inoculation  is  slowed 
(22).  Hence,  these  octarepeats  are  not  necessary  for  disease, 
but  i),  they  can  impart  structure  in  the  presence  of  copper  that 
confers  proteinase  resistance  to  the  prion  protein,  and  ii), 
their  presence  enhances  disease  incubation.  Indeed,  it  is 
known  that  the  octarepeat  region  adopts  a  unique  structural 
motif  when  copper  is  bound  (23).  We  can  rationalize  these 
data  by  assuming  the  following  affinity  ranking:  cuprizone  > 
nonoctarepeat  site  >  D-PEN  >  octarepeat  site.  This  ranking 
is  plausible  given  the  hexadentate  chelation  of  cuprizone 
versus  at  most  tetradentate  coordination  of  D-PEN.  With  this 
assumption,  D-PEN  lifts  the  proteinase  resistant  structure 
associated  with  copper  binding  to  octarepeats,  and  this  slows 
conversion  to  the  infectious  form.  These  assumptions  may  be 
tested  by  performing  competitive  binding  experiments  in 
vitro  with  D-PEN  (18),  cuprizone  (19),  and  the  relevant 
copper  binding  segments  of  the  prion  protein.  We  note  that 
the  high  affinity  of  cuprizone  for  divalent  copper  apparently 
leads  to  stripping  of  copper  from  other  proteins  so  that  the 
in  vivo  phenomena  need  not  be  limited  to  prions;  if  our  as¬ 
sumptions  are  correct,  these  nonprion  related  effects  must 
lead  to  more  rapid  degeneration  of  neurons  than  conversion 
to  PrPSc. 

To  further  test  the  mechanism,  it  is  desirable  to  perform 
mutation  studies  on  the  stretch  GGSH  or  GGTH.  The  key 
idea  is  to  mutate  out  the  H  for,  e.g.,  A,Y,G  or  other  amino 
acids  which  should  severely  mitigate  copper  binding.  We 
propose  first  that  recombinant  mutant  PrPc  from  residues 
90-230  (with,  e.g.,  H96A  in  the  human  form)  be  allowed 
to  aggregate  in  copper-full  and  copper-free  environments. 
Using  the  reduced  length  prion  will  avoid  the  confounding 
effects  from  the  octarepeat  regions  discussed  above.  If  the 
mechanism  is  correct,  the  fibril  formation  rate  should  be 
approximately  the  same  as  wild-type  fibrils  grown  in  copper- 
free  solutions.  Second,  one  can  transfect  cultured  neurons 
with  the  mutated  DNA  construct  and  look  for  increased 
susceptibility  to  infection.  Third,  assuming  positive  results 
from  the  fibril  growth  and  cultured  neuron  experiments, 
transgenic  mice  with  the  H95  mutated  away;  these  mice 
should  prove  to  have  shorter  incubation  times  for  a  given 
dose  than  wild-type  mice.  Finally,  fibrillization  experiments 
on  recombinant  marsupial  PrPc  with  and  without  copper 
should  be  informative  since  the  H  is  replaced  with  Y24;  our 
prediction  based  upon  the  present  hypothesis  would  be  that 
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the  copper  should  have  reduced  inhibition  of  fibril  formation 
relative  to  growth  of  recombinant  human  prion  protein. 
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Mammalian  prion  proteins  (PrP)  are  of  significant  public  health  interest.  Yeasts  have  proteins,  which  can 
undergo  similar  reconformation  and  aggregation  processes  to  PrP,  without  posing  a  threat  to  the  organism. 
These  yeast  “prions,”  such  as  SUP35,  are  simpler  to  experimentally  study  and  model.  Recent  in  vitro  studies 
of  the  SUP35  protein  found  long  aggregates,  pure  exponential  growth  of  the  misfolded  form,  and  a  lag  time 
which  depended  weakly  on  the  monomer  concentration.  To  explain  this  data,  we  have  extended  a  previous 
model  of  aggregation  kinetics  along  with  a  stochastic  approach.  We  assume  reconformation  only  upon  aggre¬ 
gation  and  include  aggregate  fissioning  and  an  initial  nucleation  barrier.  We  find  that  for  sufficiently  small 
nucleation  rates  or  seeding  by  a  small  number  of  preformed  nuclei,  the  models  achieve  the  requisite  exponen¬ 
tial  growth,  long  aggregates,  and  a  lag  time  which  depends  weakly  on  monomer  concentration.  The  spread  in 
aggregate  sizes  is  well  described  by  the  Weibull  distribution.  All  these  properties  point  to  the  preeminent  role 
of  fissioning  in  the  growth  of  misfolded  proteins. 

DOI:  10.1 103/PhysRevE. 72. 051915  PACS  number(s):  82.39.Fk 


I.  INTRODUCTION 

Bovine  spongiform  encephalopathy  (BSE)  in  cows, 
Scrapies  in  sheep,  Creutzfeldt-Jakob  Disease  (CJD)  and 
Kuru  in  humans  are  all  diseases  caused  by  a  specific  mis¬ 
folded  protein  residing  on  neurons  [1,2].  This  prion  protein 
(denoted  PrP'  in  its  normal  form  and  PiP"  in  its  misfolded 
form)  is  present  in  all  mammals  and  its  full  function,  to  date, 
is  still  unknown.  A  growing  body  of  evidence  strongly  im¬ 
plies  that  this  disease  propagates  not  by  nucleic  acids,  such 
as  DNA  or  RNA,  but  by  misfolded  proteins  (PrPvr)  [3,4]. 
In  our  current  understanding  of  prion  diseases,  the  autocata- 
lytic  misfolding  of  the  prion  proteins  plays  a  central  role. 
The  misfolded  form  PrP"  entices  normal  versions  of  the  cel¬ 
lular  prion  protein  (PrP'  )  to  change  conformation  to  the  mis¬ 
folded  or  disease  causing  form.  Although  PrPe  and  PrPIC 
have  the  same  amino  acid  sequence  (230  amino  acids  or 
residues),  PrPsc  has  a  higher  beta  sheet  content  than  its  nor¬ 
mal  form  [1].  In  humans,  the  disease  is  mostly  sporadic, 
perhaps  caused  by  a  rare  spontaneous  misfold  of  the  protein. 
But,  infectious  forms  are  also  known  to  occur,  such  as  in  the 
case  of  new-variant  CJD  caused  by  eating  BSE  infected 
meat.  In  the  latter  case,  infectious  agents  are  presumably 
misfolded  “seeds”  that  have  entered  the  body  from  outside. 
Developing  an  understanding  of  the  misfolding  process 
should  provide  insights  to  prevention  and/or  cure  of  these 
diseases. 

A  class  of  proteins  in  yeast  (SUP35)  undergo  a  conforma¬ 
tion  change  similar  to  the  mammalian  PrP.  However,  in  case 
of  yeast,  this  does  not  kill  the  organism.  In  fact,  reconfigured, 
aggregated  forms  lead  to  a  new  phenotype  [5].  Because  they 
are  not  toxic,  the  misfolding  is  much  easier  to  study  for  these 
proteins.  Besides,  it  has  been  shown  that  aggregates  pro¬ 
duced  in  vitro  can  lead  to  the  same  phenotype,  when  added 
to  a  yeast  cell  [6-8].  Thus,  it  is  possible  to  study  the  molecu¬ 
lar  aggregation  processes  for  yeast  prions  in  a  test  tube  in 
tremendous  detail,  and  a  large  quantity  of  experimental  data 
is  indeed  available  [9,10]. 


In  this  paper,  we  explore,  using  conventional  kinetic  as 
well  as  stochastic  models,  the  aggregation  dynamics  for 
model  yeast  prions  with  a  particular  goal  of  explaining 
in  vitro  data  [9,10]  including  (i)  linear  aggregates  which 
grow  by  monomer  addition,  (ii)  a  modest  concentration  de¬ 
pendence  of  initial  aggregate  growth  time  varying  roughly 
inversely  as  the  square  root  of  the  monomer  concentration, 
(iii)  an  observed  distribution  of  aggregate  sizes  with  a  mode 
that  is  large  (lOs-lOOs  of  monomers),  and  (iv)  a  sigmoidal 
growth  that  is  nearly  exponential.  We  extend  the  nucleation 
and  fissioning  model  of  Collins  et  al.  [9,10]  from  a 
moments-only  model  to  monitor  individual  polymer  lengths 
and  show  that  we  can  readily  obtain  the  above  listed  features 
in  the  data.  In  particular,  we  argue  against  a  model  based 
upon  rapid  equilibration  with  micelles,  which  can  also  pro¬ 
vide  a  weak  concentration  dependence  to  initial  aggregation 
times  [11,12],  but  not  the  exponential  growth  implied  by 
including  explicit  fission  of  aggregates.  The  fission  rate  plays 
a  central  role  in  determining  the  time  scales  as  well  as  the 
mean  aggregate  length.  The  range  of  aggregate  sizes  assumes 
a  Weibull  distribution,  also  well  known  from  other  processes 
where  fissioning  is  important.  Thus  the  agreement  with  ex¬ 
periments  highlights  the  importance  of  fissioning  in  the 
growth  dynamics. 

We  study  our  kinetics  both  with  continuous  time  model¬ 
ing  and  via  the  discrete  time  stochastic  Gillespie  algorithm. 
The  motivation  for  applying  the  latter  is  to  determine  where 
stochastic  effects  from  the  relatively  rare  processes  of  large 
length  polymer  generation  may  be  observable  in  experiment. 

In  considering  the  relevance  of  these  studies  to  mamma¬ 
lian  prions,  there  is  an  important  caveat  that  arises  from  the 
fact  that  mammalian  prions  are  glycosyl  phosphatidyl- 
inositol  (GPI)  anchored  on  neuron  cells  and,  thus,  move  in  a 
two-dimensional  plane  and,  hence,  can  lead  to  areal  aggre¬ 
gates  [13,14].  The  aggregation  and  fissioning  processes  are 
likely  to  be  very  different  there.  Nevertheless,  the  dominance 
of  fissioning  can  lead  to  long-time  growth  kinetics,  which  are 
in  many  ways  closely  related. 
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FIG.  1.  Aggregation  model  where  circles  are  normal  proteins 
(psi-)  and  squares  are  the  misfolded  type  [PSI+], 

II.  MODELS 

We  consider  a  one-dimensional  aggregation  model  that 
includes  a  nucleation  barrier,  elongation  by  monomer  addi¬ 
tion,  and  bssion  (Fig.  1).  The  same  kind  of  kinetics  was 
considered  by  Collins  et  al.  in  the  form  of  moments  of  the 
full  polymer  kinetic  equations,  but  not  in  the  explicit  length 
resolved  form  considered  here  [9,10].  The  assumption  that 
monomer  (rather  than  oligomer)  addition  dominates  and  the 
need  for  fission  are  underscored  by  the  data  we  seek  to  de¬ 
scribe  on  yeast  prions  [9,10].  We  shall  distinguish  the  model 
quantitatively  from  a  nonfission  based  model  later  in  the 
paper. 

A  nucleation  event  occurs  very  rarely  and  is  composed  of 
two  normal  [psi-]  proteins  interacting  in  such  a  way  that  a 
misfold  occurs  and  they  bind  together;  this  two  prion  aggre¬ 
gate  is  called  a  [PSI+]  dimer.  The  actual  size  of  the  nucleus  is 
not  a  critical  parameter  in  our  model,  as  long  as  the  nucle¬ 
ation  process  is  much  slower  than  aggregation.  Once  an  ag¬ 
gregate  has  been  created,  it  can  elongate  by  monomer  addi¬ 
tion.  Our  model  does  not  consider  fusion  between  two 
oligomers.  The  rate  of  fissioning  becomes  important  only 
when  an  oligomer  of  large  enough  size  arises.  After  an  oli¬ 
gomer  has  fissioned  into  smaller  pieces,  the  individual  pieces 
can  grow  and  fission.  When  a  fission  leads  to  a  monomeric 
product,  we  assume  that  it  returns  to  the  properly  folded 
state.  This  process  leads  to  a  steady  state  distribution  of 
oligomer  sizes  and  a  pure  exponential  growth.  It  continues 
until  the  monomer  concentration  begins  to  dwindle 
significantly. 

We  note  that  the  lag  time  for  the  in  vitro  growth  of  yeast 
prion-like  proteins  is  reproducible  in  experiments,  and  shows 
a  weak  dependence  on  monomer  concentration.  Hence,  it  is 
not  likely  that  it  can  be  associated  with  pure  nucleation  (at 
least  at  low  initial  monomer  concentration),  which  being  a 
rare  event  will  lead  to  a  distribution  of  time  scales  as  well  as 
a  rather  strong  dependence  on  monomer  concentration,  de¬ 
pending  on  the  size  of  the  nucleus.  We  will  assume  that  the 
lag  time  arises  from  the  growth  of  misfolded  material  to  an 
observability  threshold.  This  assumption  is  valid  if  the  nucle¬ 
ation  time  scales,  which  are  long  compared  to  other  micro¬ 


scopic  time  scales,  are  still  shorter  than  this  latter  time  scale. 
It  is  also  valid  when  a  small  number  of  nucleated  seeds  exist 
at  the  time  of  the  preparation  of  the  experimental  samples. 
As  shown  elsewhere,  a  weak  concentration  dependence  to 
lag  time  can  derive  from  rapid  equilibration  of  monomer 
concentration  with  nonfibrillar  oligomers  [11,12];  this  model 
is  not  supported  in  the  current  context,  and  we  shall  discuss 
this  later  in  this  section. 


A.  Kinetic  model 

Our  first  approach  is  to  use  kinetic  equations  to  describe 
the  aggregation  process  [15,16].  In  order  to  use  the  rate  equa¬ 
tions,  one  must  assume  a  large  number  of  monomers.  The 
kinetic  rate  equations  for  our  model  are  as  follows 


dnx 

dt 


=  -  2  p, 


nucn  1 


N-\  N-\ 

2p+«l2  ni  +  2pf^j  nh 
i=2  i= 2 


(1) 


dn2  2 
~  Pnuc11 1 


N- 1 

2p+n  |  n2  -  pfi2  +  2pf^  nh 
i'=3 


(2) 


dnk 

dt 


N- 1 

=  2p+nk_tn]  -  2 p+nknx  -{k-  1  )pfik  +  2pf  2)  «„ 

i-k+2 


(3) 


dnN 

—  =2p+nN_]ni  - pf(N  —  \)nN_u  (4) 

where  nx  is  the  monomer  concentration,  nk  is  the  concentra¬ 
tion  of  A'-mers,  and  N  is  the  longest  oligomer  kept  in  the 
calculation.  The  parameters  pnuc,  p+,  and  pj-  can  be  adjusted 
to  change  the  rate  of  nucleation,  fusion,  and  fission,  respec¬ 
tively.  Since  one  of  these  parameters  can  be  absorbed  into 
the  definition  of  the  time  f,  we  will  set  p+=  1.  Furthermore, 
we  will  set  the  initial  monomer  concentration  to  unity.  The 
equations  listed  above  form  a  system  of  A'-coupled  differen¬ 
tial  equations  which  cannot  be  solved  analytically  because  of 
the  nonlinear  terms. 

An  analytical  solution  is  possible  if  nx  can  be  assumed  to 
be  a  constant  and  n2  and  higher  are  much,  much  smaller  than 
the  monomer  concentration.  Then  the  set  of  equations  be¬ 
come  linear  and  can  be  solved  by  several  techniques,  such  as 
a  Laplace  transformation.  This  is  not  a  bad  approximation, 
for,  as  we  will  see,  during  much  of  the  growth  process  the 
monomer  concentration  is  nearly  constant.  Only  at  the  end, 
the  monomer  concentration  begins  to  dwindle  and  the 
amount  of  misfolded  monomers  saturates.  It  is  useful  to  de¬ 
fine  the  zeroth  and  first  moments  of  the  aggregate  size  dis¬ 
tribution  as 


N 

A  =  2  nh 

i= 2 


(5) 
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N 

M  =  2  *  *  nn  (6) 

1=2 

where  A  is  the  number  of  polymer  aggregates  and  M  is  the 
total  number  of  monomers  in  aggregate.  Since  the  overall 
number  of  proteins  is  conserved,  the  rate  of  change  of  M  is 
given  by 

dM  dn  | 

- =  -— 7 

dt  dt 

whereas  the  rate  of  change  of  A  is  given  (ignoring  the  cutoff 
N)  by  the  equation 

dA  -)  ,  ■ 

—  =Pnucn\  -Pfll+Pflj  (;-3)«/.  (8) 

i= 3 


If  we  ignore  the  terms  depending  on  n2  and  assume  « ,  is  a 
constant,  we  get  two  linear  coupled  equations, 

dA 

—  =  const  +  pfM  -3pfA,  (9) 


dM 

- =  const  +  (2  p+nx 

dt 


-2 Pf)A. 


(10) 


We  set  p+=  1  to  set  the  unit  of  time.  Anticipating  M  > A , 
which  means  that  the  mean  aggregate  size  is  much  bigger 
than  unity.  The  above  equations  mean  that  the  M  and  A  both 
grow  exponentially  as  ekl,  with  X  given  by 

X  =  V 2pf{nx  -pf)~  i2pfn x .  (11) 

During  the  exponential  growth,  the  mean  aggregate  size  is 
given  by 


If  we  further  assume  that  the  lag  time  is  given  by  the  time  it 
takes  for  the  misfolded  proteins  to  reach  a  detectable  thresh¬ 
old  Mj-,  then  the  lag  time  is  given  by 


In 


Mr 


‘lag  ' 


C(n  i)2 
V2  pfnl 


(13) 


where  C  is  a  constant  and  the  n  x  dependence  inside  the  loga¬ 
rithm  comes  from  the  prefactors.  We  will  see  that  these  re¬ 
sults  are  confirmed  by  a  complete  numerical  integration  of 
the  differential  equations.  Thus,  this  model  demonstrates 
the  requisite  /-lag  dependence  seen  in  the  experimental 
data,  namely  the  inverse  square  root  dependence  on  the 
monomer  concentration,  up  to  logarithms,  shown  in  Fig.  2. 
We  note  that  the  fit  in  Fig.  2  is  not  sensitive  to  the  logarith¬ 
mic  terms. 


Kinetic  model  results 

A  numerical  integration  approach  was  used  to  study  the 
full  set  of  equations.  We  used  the  fourth  order  Runge-Kutta 
algorithm  for  solving  differential  equations  and  to  obtain 


FIG.  2.  Lag  time  ( ttag )  vs  initial  monomer  concentration  [/z x (0)]. 
The  crosses  are  data  from  Collins  et  al.  in  Ref.  [9].  The  continous 
curve  is  a  fit  to  the  equation  f/a£=«o[oi(0)]_I/2  log  axhix(t)),  with 
ci0=9.22  min* (pM)m  and  fl[  =  1.80X  10 6pM,  following  Eq.  (13). 

length  distributions  and  growth  curves  for  the  system. 

In  Fig.  3,  we  show  plots  of  log(M)  as  a  function  of  time. 
We  have  chosen  a  set  of  parameters  yielding  mean  lengths  in 
the  10.?- 100s  and  displaying  manifest  exponential  growth  in 
M{t)  vs  t  plots.  It  is  evident  that  there  is  a  regime  where 
log (M)  varies  linearly  with  time,  implying  a  pure  exponen¬ 
tial  growth.  This  region  of  exponential  growth  is  not  limited 
to  just  a  few  chosen  parameter  values,  but  can  be  achieved 
over  a  wide  range  of  pf  values.  In  Fig.  3,  we  plot  the  total 
number  in  aggregate  vs  time,  now  on  a  linear  scale.  The 
sigmoidal  growth  is  now  evident  as  the  exponential  increase 
tapers  off  when  the  monomers  begin  to  deplete.  This  is  also 
consistent  with  the  experimental  results  [9]. 

As  discussed  earlier,  the  mean  length  of  the  aggregates 
has  a  simple  inverse  square-root  dependence  on  the  rate  of 
fission.  Thus  by  tuning  the  rate  of  fission,  it  is  possible  to 
obtain  mean  aggregate  lengths  (see  Fig.  4),  which  appear 


t>0 

O 


FIG.  3.  Log[M(t)]  vs  time,  with  parameter  values  n1(0)  =  l,  p+ 
=  1,  p,mc=  10-8.  The  linear  regions  illustrate  regions  of  pure  expo¬ 
nential  growth.  M(t)  vs  time  illustrating  sigmoidal  growth  is  shown 
in  the  inset.  Note  that  time  is  in  units  of  1  /[/i](0)p+]. 


051915-3 


26 


KUNES,  COX,  AND  SINGH 


PHYSICAL  REVIEW  E  72,  051915  (2005) 
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FIG.  4.  Length  distribution,  nt  vs  length  L, 
from  the  kinetic  equations  during  the  exponential 
growth  phase,  compared  with  the  Weibull  distri¬ 
bution,  Eq.  (14). 


qualitatively  consistent  with  those  observed  experimentally 
displaying  lengths  in  the  tens  of  hundreds  (see  especially 
slide  9  at  Ref.  [10]).  In  Fig.  5,  this  mean  length  is  plotted  vs 
time  for  a  different  fission  rate  pf.  We  see  that  the  mean 
aggregate  length  is  constant  during  the  exponential  growth. 
In  Fig.  6,  we  show  the  variation  of  the  exponential  growth 
rate  with  compared  with  the  analytic  calculations.  One  can 
see  that  they  agree  closely.  Note  also  that  pf  values  of  order 
10~4  will  give  the  lengths  seen  in  vitro. 

We  can  also  determine  the  distribution  of  aggregate  sizes 
in  our  calculations.  The  concentration  of  aggregates  of  dif¬ 
ferent  sizes  relative  to  the  initial  concentration  of  monomers 
are  shown  in  Fig.  4  for  a  given  value  of  the  fission  parameter 
(j>f=  4X  1(T4)  at  three  different  times  during  the  exponential 
growth  phase  of  the  simulation.  Clearly,  the  distribution  is 
fairly  stable  during  the  exponential  growth  phase,  further 
strengthening  the  argument  for  a  simple  steady  state.  These 
length  distributions  are  well  fitted  by  a  Weibull  distribution 
given  by  the  relation. 


FIG.  5.  Mean  length  vs  time  for  different  pf,  «j(0)=  1,  p+=  1, 
Pnuc=  10-8.  The  plateau  corresponds  to  the  constant  mean  length 
during  the  exponential  growth  phase.  Note  that  time  is  in  units  of 

l/[«i(0)p+]. 


fix)  = 


Oq 
Cl  i 


(14) 


Apart  from  an  overall  normalization,  the  two  key  parameters 
of  the  distribution  are  the  mean  ai  and  the  power  law  at 
small  x  set  by  ciq  -  1.  This  distribution  is  found  to  be  ubiqui¬ 
tous  in  nature  [17].  For  example,  a  dropped  piece  of  coal  will 
have  shattered  pieces  that  follow  an  asymmetrical  distribu¬ 
tion  with  more  pieces  on  the  smaller  end.  Our  model  follows 
a  similar  idea  (i.e.,  taking  a  larger  length  and  shattering  it 
into  smaller  lengths)  and  a  similar  distribution.  In  our  simu¬ 
lations,  the  quantity  a0  is  close  to  2  for  a  range  of  parameters 
studied.  During  subsequent  times,  the  peak  of  this  curve 
shifts  to  the  left  as  saturation  occurs  and  the  only  process  left 
is  to  fission.  (The  shift  in  peak  at  later  times  are  not  shown 
here.) 


FIG.  6.  Exponential  growth  rate  \  vs  pt  compared  with  simple 
square-root  dependence  predicted  by  Eq.  (11)  for  different  param¬ 
eter  values,  nj(0)  =  l,  p+=  1,  pnuc=  10  . 
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B.  Comparison  to  nonfissioning  micelle  model 

The  coupling  of  the  observation  of  weak  concentration 
dependence  with  exponential  growth  strongly  supports  a 
model  dominated  by  fissioning  of  aggregates  constructed  by 
monomer  addition.  Indeed,  centrifugation  studies  confirm 
that  monomers  are  the  dominant  species  in  the  in  vitro  me¬ 
dium,  and  examination  of  the  aggregation  kinetics  suggests  a 
critical  nucleus  of  small  size  [9], 

It  is  thus  important  to  contrast  with  other  models  for 
which  a  weak  concentration  dependence  can  be  derived.  No¬ 
tably,  the  work  of  Lomakin  et  al.  [11,12]  considers  kinetics 
with  monomer  growth  of  polymers  as  well  as  micellar  oligo¬ 
mers  composed  of  A4  monomers  of  concentration  nM  which 
may  or  may  not  be  off  pathway.  This  model  may  be  of  rel¬ 
evance  to  other  in  vitro  studies  of  yeast  prion  aggregation 
[18]  for  which  there  is  evidence  of  oligomer  formation  and 
growth  and  a  weak  concentration  dependence  to  the  lag 
time.  If  the  oligomers  are  off  pathway,  then  assuming  an 
equilibration  with  monomers  leads  to  a  buffering  of  the 
monomer  concentration  assuming  an  initially  high  input  con¬ 
centration  and  sufficiently  rapid  equilibration.  There  is  a 
crossover  or  critical  concentration  n*  for  this  dictated  by  the 
equilibrium  constant  KM  between  monomers  and  micelles 
(■ KM=nM/n given  by 

n=(MKM)m-M.  (15) 

Lomakin  et  al.  argue  that  for  a  total  equilibrated  protein  con¬ 
centration  nT  obtained  prior  to  fibril  growth, 


=  nl  +  MnM 

(16) 

-•cr 

(17) 

From  these  equations,  it  is  clear  that  there  is  an  approximate 
two-state  behavior  provided  J\Ai>  1.  For  nT<n*,  then  nT 
—  Hi  and  aggregation  kinetics  can  be  strongly  concentration 
dependent.  However,  for  nT^n*,  then  n ,  ~  n* ,  that  is,  the 
micelles  buffer  the  monomer  concentration.  This  would  lead 
to  a  weak  dependence  of  fibrillar  growth  upon  subsequent 
increase  of  initial  monomer  concentrations. 

There  are  two  routes  by  which  this  model  can  be  applied 
to  the  data  of  Collins  et  al.  [9,10].  First,  if  the  fibrils  are 
formed  from  monomer  addition  to  the  micelles  without  re¬ 
sulting  fission  of  fibrils,  then  one  will  expect  a  quadratic 
growth  in  time  [19,20]  to  the  total  number  of  aggregated 
proteins,  namely 


M{t)~-p2+n\nMt2,  (18) 

~  ^^pln]+M(n*)'~MnMt2,  (19) 

which  is  very  weakly  dependent  upon  «i(0)  provided  «i(0) 
~nT^n*.  but  varies  strongly  as  n1(0)2+A(  in  the  opposite 
limit.  Solving  for  where  M(tlay)  =  M c  characteristic  of  the  ob¬ 
servability  threshold  yields 


i 

,ag  p+nx  V  (m)M 


(20) 


For  ii\{ 0)<§«*,  t!ag  varies  as  nl(0)^1+M,2\  at  odds  with  ex¬ 
periment  [9],  while  for  «i(0)  and  A4  >  1 ,  tIag  is  approxi¬ 
mately  independent  of  /q(0)  [it  will  decrease  weakly  with 
increasing  «i(0)]. 

Second,  if  the  micelles  are  off  pathway,  then  their  princi¬ 
pal  impact  on  the  fibrillar  aggregation  is  on  buffering.  As¬ 
suming  a  critical  fibrillar  nucleus  of  length  p,  which  equili¬ 
brates  with  the  monomers  via  a  constant  Kp=nplnpx ,  the  same 
considerations  of  aggregation  kinetics  as  in  the  preceding 
paragraph  give  a  lag  time  given  by 


1  /  2  Mc 

P+n\  V  Kp{nx)p 


(21) 


For  n1(0)<^n*,  this  varies  as  n,(0)  ,1+A2)  js  at  odds  with 
Collins  et  al.  [9],  while  it  is  approximately  independent  of 
n^O)  for  /i](0) and  p> 2.  For  small  nuclei  (p  —  2-5, 
say),  the  dependence  can  be  relatively  strong  on  7q(0)  even 
in  this  limit. 

Hence,  via  either  on-pathway  or  off-pathway  micelles,  in 
this  model  there  is  a  route  to  weak  dependence  of  lag  time 
(defined  as  the  threshold  for  observation)  upon  initial  mono¬ 
mer  concentration  requiring  nj(0)  and  A4  ,p  >  1 .  How¬ 
ever,  there  are  several  problems  in  connecting  this  to  the 
experiments  of  Collins  et  al.  First,  exponential  growth  is 
unambiguously  observed,  indicating  that  fission  of  fibrils  oc¬ 
curs.  The  micelle  model  of  this  subsection  has  no  fission. 
This  highlights  the  important  role  of  fission  for  the  prion 
phenomenon.  Second,  in  these  experiments,  the  starting  so¬ 
lution  contains  almost  entirely  monomers  and  not  micelles. 
This  implies  that  even  if  micelles  form,  one  is  always  in  the 
limit  of  ni(0)<n",  leading  to  strong  concentration  depen¬ 
dence.  Third,  the  width  of  the  fibrils  obtained  in  this  experi¬ 
ment  is  apparently  monomeric,  which  is  not  readily  compat¬ 
ible  with  growth  by  monomer  addition  to  micelles,  and 
rather  must  be  obtained  by  addition  to  a  critical  fibril 
nucleus.  The  data  suggest  that  this  critical  fibril  nucleus  is  of 
length  p^6,  not  consistent  with  the  assumption  p  >  \ . 

Hence,  we  do  not  believe  that  the  micelle  based  model 
can  explain  the  data  of  Collins  et  al.  [9,10]. 


C.  Stochastic  model 

In  this  section,  we  develop  a  stochastic  treatment  of  the 
model,  similar  to  a  model  proposed  by  Poschel  et  al.  [21]. 
Our  main  motivation  is  to  be  able  to  treat  species  with  small 
numbers,  in  which  case  the  continuous  deterministic  ap¬ 
proach  will  break  down.  As  we  will  see,  the  stochastic  ap¬ 
proach  largely  agrees  with  the  deterministic  model  for  the 
parameter  values  considered  here  for  aggregate  quantities 
[e.g.,  M(t)].  However,  for  distributions  there  can  be  substan¬ 
tial  stochastic  noise  at  short  times.  Another  advantage  of  the 
stochastic  approach  is  that  it  can  be  readily  extended  to  study 
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TABLE  I.  Propensities  of  stochastic  processes. 


Fusion  propensity  for 

monomer  addition 

Fission  propensity  for  monomer 

Higher  order 
fission  term 

propensity 

Chain  length  2 

Chain  length  3 

fusion  term 

Chain  length  2 

Chain  length  3 

gnucNMl-l) 

2gJfiN2 

2g+NjN3 

gfN2 

2g/N3 

two-dimensional  aggregates  as  well  as  the  problem  of  mul¬ 
tiple  prion  “strains,”  in  which  case  the  treatment  of  a  small 
number  of  heteroaggregates  would  be  important  to  model  the 
extent  to  which  strains  breed  true. 

The  Gillespie  algorithm  provides  an  exact  way  to  treat  the 
stochastic  problem  of  chemical  reactions  [22].  In  our  case, 
the  polymers  of  different  lengths  are  the  different  chemical 
species.  The  processes  of  nucleation,  monomer  addition,  and 
fission  are  assumed  to  be  stochastic.  In  other  words,  the 
number  of  polymers  of  different  length  at  time  t,  only  define 
the  propensities  (or  normalized  probability)  for  the  different 
reactions  to  happen  at  that  time.  Once  a  reaction  takes  place, 
the  number  of  polymers  is  altered  and  the  propensities  are 
changed.  The  Gillespie  algorithm  is  a  Monte  Carlo  treatment 
that  deals  with  the  stochastic  process  by  using  a  pair  of  ran¬ 
dom  numbers  at  each  step,  one  to  decide  which  event  will 
occur  next  and  another  to  decide  how  long  will  it  be  until  the 
next  event  takes  place.  This  process  can  be  repeated  to  fol¬ 
low  the  dynamical  behavior  of  the  system. 

Our  model  has  the  corresponding  propensities 


gnuMNx  -  1)  +  2g+N Ni  +  gjZ  (i  -  l)Nj,  (22) 


i=2 


where  gnuc,  g+,  and  gf  are,  respectively,  the  nucleation,  fu¬ 
sion,  and  fission  parameters  in  the  stochastic  model.  In  our 
model,  each  propensity  was  assigned  a  bin  in  an  array.  A 
random  number  generator  decides  which  bin  is  selected. 

Table  I  illustrates  physically  the  propensities  at  a  given 
time.  A  normalized  random  number  generator  selects  which 
element  in  the  array  will  occur  and  the  time  it  took  to  create 
that  event.  After  an  event  occurs,  the  corresponding  propen¬ 
sities  are  updated  and  the  process  is  repeated. 

In  order  to  compare  the  results  achieved  from  the  continu¬ 
ous  to  stochastic  model,  a  mapping  of  the  Gillespie  param¬ 
eters  to  the  kinetic  parameters  is  required.  For  this,  we  need 
to  develop  an  approximate  equation  satisfied  by  the  mono¬ 
mer  concentration  in  the  stochastic  simulations  and  compare 
it  with  the  rate  equations.  In  the  stochastic  case,  one  has  a 
master  equation  that  relates  the  probability  distribution  asso¬ 
ciated  with  a  different  number  of  polymers  at  time  t+dt  to 
those  at  time  t.  We  will  make  a  mean-field  approximation  for 
different  kind  of  k-mers,  namely  ( V; * Nk) = (A,)  *(Nk).  Thus, 
we  can  consider  the  approximate  master  equation,  which 
only  tracks  changes  in  the  monomer  number.  We  arrive  at  the 
equation. 


Pl{Nl,t  +  dt)  =  Pl{Nl,t) 


1  ~  8, mcNiin^  -  l)dt-g+N Njdt 


i=2 


■  2 gf^j  Njdt 


i=2 


+  2P,(A,  -  1  ,t)g£  Njdt 


i= 3 


+  Pt(Nt  -  2,t)gfN2dt  +  H|(A'|  +  1  .rWVj 


x2  Nidt  +  P1(Nl  +  2,t)gnucN1(Ni  -  1  )dt, 

/=2 

(23) 

where  A]  is  the  number  of  monomers,  P\(N\ , t+dt)  is  the 
probability  to  have  Nx  monomers  at  a  later  time  t+dt. 
P|(V|  ,t)  is  the  probability  to  have  N1  monomers  at  time  t.  It 
is  multiplied  by  the  probability  that  no  reaction  occurs  in 
time  dt  that  changes  A).  The  rest  of  the  terms  represent  the 
probability  to  have  a  different  number  of  monomers  at  time  t 
but  then  a  reaction  happens  in  time  dt  leading  to  Ai  mono¬ 
mers  at  time  t+dt. 

Now  using  the  definition  of  the  derivative,  one  can  rear¬ 
range  the  equation  and  multiply  both  sides  by  Aj  to  get  the 
average  rate  of  change  of  Nt .  After  shifting  some  indices  to 
get  every  probability  to  have  the  form  P\(Nj  ,t)  and  using  the 
appropriate  volume  element  (denoted  by  V)  to  normalize  the 
counts  to  a  concentration,  one  arrives  at  the  equation  for  the 
mean  number  of  monomers. 


1  dnx 
Vg+dtgil 


=  ~  2gnucn\  -  2n^  n,  +  2^2 


i=2 


V 


n-i ■ 


(24) 


By  comparing  with  the  rate  equations  discussed  before,  we 
obtain  the  parameter  mapping 


P+=  1> 

(25) 

Pnuc  ~  Snuct 

(26) 

ii 

(27) 

t  kinetic  8  +  ^ gillespie  • 

(28) 

Results  stochastic  approach 

After  achieving  the  parameter  mapping,  we  studied  the 
stochastic  models  with  106  initial  monomers.  The  numerical 
results  are  shown  in  the  next  few  figures.  In  Fig.  7,  the 
growth  of  the  aggregate  material  is  shown  on  a  logarithmic 
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FIG.  7.  Log[M(f)]  and  M(t),  shown  in  the  inset,  vs  time  ob¬ 
tained  by  the  Gillespie  algorithm  with  IVi(0)  =  10  ,  g+=l,  gnuc 
=  10“8  for  different  fission  rate  gt.  Note  that  time  is  in  units  of  1  /g+. 

and  on  a  linear  scale.  It  is  evident  that  they  are  in  excellent 
agreement  with  the  corresponding  calculations  for  the  kinetic 
model.  In  Fig.  8,  we  show  plots  of  the  mean  length  as  a 
function  of  time,  and  in  Fig.  9,  the  corresponding  distribution 
of  length  scales  during  the  exponential  growth  phase.  It  is 
again  evident  that  there  is  a  steady  state  during  the  exponen¬ 
tial  growth  and  the  mean  lengths  reach  a  plateau  value.  The 
stochastic  effects  are  much  larger  in  the  distribution  and  in 
the  mean  length,  but  the  overall  results  agree  well  with  the 
kinetic  model.  The  fits  to  the  Weibull  distribution  are  again 
good  with  comparable  parameters.  The  differences  between 
the  kinetic  and  stochastic  Weibull  parameters  ( a0  and  «,)  are 
less  than  10%.  This  shows  that  for  106  monomers  and  for 
mean  aggregate  lengths  up  to  a  few  hundred,  the  dominant 
species  still  occur  in  large  enough  numbers  so  that  stochastic 
effects  do  not  change  the  results  in  a  significant  way. 

III.  DISCUSSION 


In  this  paper,  we  have  used  a  deterministic  kinetic  and  a 
stochastic  Monte  Carlo  approach  to  model  the  dynamics  of 


FIG.  8.  Mean  length  vs  time  in  units  of  1  / g+  for  different  fission 
rates  gt  obtained  by  the  Gillespie  algorithm.  The  parameter  values 
are  A1(0)  =  106,  g+=  L  gnUc=  FT8. 


FIG.  9.  Length  distributions  at  different  times  during  the  expo¬ 
nential  growth  phase  obtained  by  the  Gillespie  algorithm  for  g, 
=  400  and  fitted  by  the  Weibull  distribution. 

yeast  prion  growth.  For  the  parameters  relevant  to  in  vitro 
experiments,  where  mean  aggregate  size  is  of  order  lOOs  of 
monomers,  and  the  number  of  monomers  is  larger  than  106, 
the  two  approaches  agree  well  with  each  other,  showing  that 
stochastic  effects  are  not  dominant.  The  studies  confirm  that 
the  key  experimental  features  of  the  growth  of  yeast  prions 
can  be  well  captured  by  these  nucleation-growth-hssioning 
models.  They  lead  to  an  exponential  and/or  sigmoidal 
growth,  an  inverse  square-root  dependence  of  the  lag  time  on 
the  monomer  concentration,  and  aggregate  sizes  that  depend 
on  the  rate  of  fission.  We  predict  that  the  aggregate  size 
distributions  should  be  Weibullian,  reflecting  the  importance 
of  fissioning  in  the  growth  process.  We  also  predict,  for  suf¬ 
ficiently  large  times,  that  the  steady  state  aggregate  length 
will  drop  as  fissioning  converts  large  aggregates  to  those  of 
minimal  size.  We  have  argued  that  models  in  which  mono¬ 
mers  equilibrate  with  micelles  but  do  not  fission  [11,12], 
which  are  capable  of  producing  lag  times  having  a  weak 
concentration  dependence  for  a  sufficiently  high  initial 
monomer  concentration,  are  in  fact  not  appropriate  for  the 
data  of  Collins  et  al.  [9,10]. 

We  note  that  the  in  vitro  aggregation  work  of  Serio  et  al. 
[18]  appears  to  come  to  different  conclusions  than  that  of 
Collins  et  al.  [9]  and  may  be  ripe  for  a  discussion  in  terms  of 
the  micelle  models.  We  do  not  understand  the  discrepancies 
between  these  two  sets  of  experiments. 

We  note  that  the  results  on  yeast  prions  discussed  here 
and  recent  work  on  mammalian  prions  stripped  of  their  GPI 
membrane  anchors  [23,24]  suggest  that  a  fundamentally  dif¬ 
ferent  mechanism  for  fission  is  at  play  for  the  latter  proteins 
vis-a-vis  yeast.  In  the  recent  experiments,  transgenic  mice 
expressing  the  cellular  prion  protein  without  the  residues 
necessary  for  the  GPI  anchor  are  exposed  to  anchorless  in¬ 
fectious  prions;  in  time  course  experiments,  prion  aggregates 
are  produced  which  retain  infectivity.  However,  the  infected 
mice  with  anchorless  PrPc  do  not  show  clinical  symptoms. 
Moreover,  the  time  course  data  [see  Fig.  1(c)  of  Chesebro 
et  al.  [24]}  show  two  remarkable  characteristics:  (i)  they  are 
nonsigmoidal  in  shape,  with  the  infectious  prion  content  at 
long  times  significantly  exceeding  that  of  infected  wild  type 


051915-7 


30 


KUNES,  COX,  AND  SINGH 


PHYSICAL  REVIEW  E  72,  051915  (2005) 


mice  at  death,  (ii)  the  time  scales  to  reach  the  levels  of  in- 
fectivity  characteristic  of  wild  type  mice  at  the  correspond¬ 
ing  infectious  dose  are  quite  long.  Indeed,  the  time  course 
data  can  be  fit  roughly  by  a  quadratic  in  time  curve  charac¬ 
teristic  of  nonfissioning  aggregation  to  an  initial  concentra¬ 
tion  of  infectious  seeds.  Evidently,  the  binding  of  cellular 
mammalian  prions  to  the  membrane  is  critical  to  the  fission¬ 
ing  process,  while  for  yeast  prions  fissioning  of  aggregates 
in  vitro  is  observed  [9]. 

In  the  future,  we  hope  to  extend  this  model  to  study  yeast 
prion  strains,  where  we  expect  the  stochastic  treatment  to  be 
the  key  to  dealing  with  rare  heteroaggregates.  An  under¬ 


standing  of  this  process  should  lead  to  insights  into  the  im¬ 
portant  problem  of  strain  dynamics  in  mammalian  prions. 
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Prion  Disease:  Exponential  Growth  Requires  Membrane  Binding 
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ABSTRACT  A  hallmark  feature  of  prions,  whether  in  mammals  or  yeast  and  fungi,  is  exponential  growth  associated  with 
fission  or  autocatalysis  of  protein  aggregates.  We  have  employed  a  rigorous  kinetic  analysis  to  recent  data  from  transgenic 
mice  lacking  a  glycosylphosphatidylinositol  membrane  anchor  to  the  normal  cellular  PrPc  protein,  which  show  that  toxicity 
requires  the  membrane  binding.  We  find  as  well  that  the  membrane  is  necessary  for  exponential  growth  of  prion  aggregates; 
without  it,  the  kinetics  is  simply  the  quadratic-in-time  growth  characteristic  of  linear  elongation  as  observed  frequently  in  in  vitro 
amyloid  growth  experiments  with  other  proteins.  This  requires  both:  i),  a  substantial  intercellular  concentration  of  anchorless 
PrPc,  and  ii),  a  concentration  of  small  scrapies  seeding  aggregates  from  the  inoculum,  which  remains  relatively  constant  with 
time  and  exceeds  the  concentration  of  large  polymeric  aggregates.  We  also  can  explain  via  this  analysis  why  mice 
heterozygous  for  the  anchor-full/anchor-free  PrPc  proteins  have  more  rapid  incubation  than  mice  heterozygous  for  anchor-full/ 
null  PrPc,  and  contrast  the  mammalian  membrane  associated  fission  or  autocatalysis  with  the  membrane  free  fission  of  yeast 
and  fungal  prions. 
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Prions  are  distinguished  from  other  amyloid  diseases  both 
by  their  infectious  character  and  the  observed  exponential 
growth  of  infectious  material  in  vivo  (1).  There  is  a  corre¬ 
spondence  to  this  of  prion-like  proteins  in  yeast  and  fungi, 
for  which  spontaneous  fission  is  reported  in  vitro  (2).  In  fact, 
it  has  been  argued  that  the  replication  necessary  for  infection 
in  mammals  and  non-Mendelian  inheritance  in  yeast/fungi 
requires  the  fission  or  autocatalysis  that  drives  the  exponen¬ 
tial  growth  (1).  It  has  been  established  for  yeast  prions  that 
additional  chaperone  proteins  most  likely  facilitate  the  fis¬ 
sion  of  aggregates  in  living  cells  ( 1 ).  It  is  an  open  question 
what  mechanism  drives  the  exponential  growth  in  mammals. 
Here  we  show  by  a  rigorous  kinetic  analysis  of  recent  disease 
time  course  data  that  the  exponential  growth  is  tied  to  mem¬ 
brane  anchoring  of  the  prion  protein,  suggesting  that  either 
mechanical  fission  of  areal  prion  aggregates  or  oligomeric 
autocatalysis  of  membrane  bound  prions  explain  the  ob¬ 
served  behavior. 

Chesebro  et  al.  (3)  recently  studied  transgenic  (Tg)  mice 
lacking  a  GPI  membrane  anchor  to  the  normal  cellular  PrPc 
protein  and  discovered  that  these  mice  grew  infectious  prions 
without  suffering  neuronal  death.  We  denote  these  anchor¬ 
less  cellular  prions  as  PrPjg  and  anchor-full  wild-type  (WT) 
cellular  prions  by  PrP^,T.  When  inoculated  with  infectious 
scrapies  prions  (PrPSc)  at  a  dose  that  induces  clinical  symp¬ 
toms  within  140-160  days  for  WT  mice,  the  Tg  mice  were 
symptomless  up  to  400-600  days,  even  though  proteinase 
resistant  PrP-res,  an  indicator  of  infectivity,  accumulated  and 
surpassed  the  maximal  WT  levels. 

In  Fig.  1  the  Tg  mice  PrP-res  concentration  ( crosses )  of 
Chesebro  et  al.  (3)  are  plotted  versus  the  square  of  time,  together 
with  a  linear  regression  fit  (line)  with  a  high  regression 
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coefficient  ( R  =  0.97).  This  time  dependence  is  consistent 
with  short  time  kinetics  described  by  linear  polymer  elon¬ 
gation  via  monomer  addition  without  fission  or  autocatalysis 
(4),  illustrated  schematically  in  Fig.  2.  Assuming  the  PrP-res 
concentration  to  be  a  proxy  for  the  total  protein  content  in 
aggregate,  simple  kinetic  arguments  predict  a  behavior,  before 
monomer  depletion  and  seed  nucleus  depletion,  of 

[PrP-res]  (t)  «  (1/2)0+  [PrP?g])2[PrPf  ]  t\  (1) 

with  p+  [PrP!f  ]  the  elongation  rate  at  intercellular  monomer 
Tg  prion  concentration  |  PrPc[  g  ] ,  and  [PrP[jc]  is  the  intercel¬ 
lular  concentration  of  seeding  nuclei  from  inoculated  scra¬ 
pies  protein  after  initial  hydrodynamic  clearance.  The 
validity  of  Eq.  1  at  long  times  suggests  that:  i),  there  is  a 
substantial  homeostatic  concentration  of  intercellular  PrP(|g 
presumably  due  to  slow  clearance,  and  ii),  |  PrP[]c]  is  hardly 
changed  implying  either  that  only  a  small  fraction  of  seeds 
grow  into  large  polymers  or  [PrP^c]  is  maintained  by  steady 
proteolytic  degradation  of  large  remnant  aggregates  from  the 
dose.  Given  a  similar  de  novo  production  rate  of  PrP^,T  in 
inoculated  WT  mice,  we  speculate  that  the  associated  sat¬ 
uration  of  [PrP-res]  arises  from  loss  of  PrP(NV  I  after  cell  death. 

This  elongation  hypothesis  is  testable  by:  i),  genetically 
engineering  mice  to  overexpress  PrP(|  g ,  which  will  quadrat- 
ically  modulate  the  PrP-res  concentration  (4),  and  ii),  by 
varying  the  initial  dose  of  PrPSc,  which  will  linearly  mod¬ 
ulate  the  PrP-res  concentration. 
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FIGURE  1  Quadratic-in-time  fit  to  infectious  prion  time  course 
data  of  Chesebro  et  al.  (3). 

Another  striking  observation  of  Chesebro  et  al.  (3)  was 
that  mice  heterozygous  for  expression  of  PrP^g  and  PrP^v| 
display  shorter  incubation  time  upon  inoculation  than  mice 
with  one  PrP(NVT  copy  and  one  inactive  fusion  construct.  We 
argue  that  this  is  due  to  an  enhanced  concentration  for  PrP(l  g 
relative  to  PrP(NV  |  and  that  PrP-res  obtained  from  PrP^g  also 
templates  PrP^j.  conversion.  The  latter  is  speculative  but 
should  be  testable. 

To  establish  the  expectation  that  the  PrP^,T  concentration 
is  lower  than  PrP^g  in  the  WT/Tg  heterozygotes,  it  is  suf¬ 
ficient  to  establish  that  PrP(NV  I  concentration  in  WT  mice  is 
lower  than  PrP.],  concentration  in  the  anchorless  Tg  mice.  To 
make  this  clear,  we  compare  estimates  for  cellular  prion 
concentrations  in  homozygous  WT  mice  with  homozygous 
Tg  mice.  Before  the  postulated  cell-death  driven  saturation 
of  infectious  material,  WT  mice  inoculated  with  a  concen¬ 
tration  [PrP^c]  of  scrapies  seeds  will  have  a  time  course  of  the 
general  form  (4) 

[PrP-res]  (f)  °c  A  +  (cosh[Ar]  —  1 ) ,  (2) 

where  A  =  ln(2)/f2  is  the  percentage  growth  rate  and  t2  is  the 
doubling  time.  The  coefficient  A+  can  be  determined  from 
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FIGURE  2  Schematic  model  for  linear  elongation  driven  growth 
of  infectious  prion  material  from  inoculated  seeds  and  anchor¬ 
less  cellular  prion  proteins. 


the  short  time  behavior  of  Eq.  2,  which,  like  Eq.  1,  is 
described  by  linear  elongation  given  by 

1/2A  +  (At)2  =  l/2(p+  [PrPwT])2[PrPnCk2,  (3) 

where  [PrP^v|  ]  is  the  homeostatic  concentration  of  mem¬ 
brane  bound  WT  PrPc  and  we  have  assumed  that  the  WT  and 
Tg  mice  have  the  same  PrP-res  elongation  coefficient  p+. 

On  the  other  hand,  at  long  times  (but  before  saturation) 
Eq.  2  gives 


[PrP-res]  (t)  x  1/2A+  exp(Ar).  (4) 

Now,  we  set  t+  =  A'd'  f2  as  the  time  it  takes  WT  mice  to 
reach  clinically  detectable  levels  of  PrP-res  concentration  at 
the  inoculum  level  generating  a  seed  nuclei  concentration 
[PrpSc],  wjlere  /\r  is  the  number  of  doublings  experienced  in 
that  process,  and  t_  is  the  time  it  takes  Tg  mice  to  reach  the 
same  clinical  concentration  of  PrP-res  for  the  same  initial  in¬ 
oculation  dose.  By  taking  suitable  ratios  to  eliminate  A+, 
[PrP®c],  and  p+,  the  ratio  of  homeostatic  concentrations  of 
cellular  prions  from  the  WT  mice  to  the  Tg  mice  is  given  by 

[PrPw-r]  /  [PrPig]  =  (l/2)ln(2)  Nd+  (f_/f  + ) 
Xexp(-ln(2)(/Vd+  /2)). 

From  Chesebro  et  al.  (3),  t+  =  150  days,  and  t_  =  400 
days.  A  reasonable  estimate  (5)  for  the  number  of  doublings 
is  N A  =  20  for  the  dose  of  Chien  et  al.  (1).  With  these 
numbers,  Eq.  5  gives  a  concentration  ratio  (and  hence 
elongation  rate  ratio)  for  WT/Tg  mice  of  0.036.  This  is 
reasonable  given  that  likely  slower  PrPc  clearance  in  the  Tg 
case  will  lead  to  a  higher  extracellular  concentration  of 
cellular  prion  protein.  By  employing  the  arguments  of  Chien 
et  al.  (1)  we  obtain  elongation  rate  values  of  0.13/day(Wt) 
and  3.5/day  (Tg).  The  former  is  in  good  agreement  with  es¬ 
timates  made  elsewhere  for  linear  elongation  based  upon 
analysis  of  dose-incubation  curves  (6). 

Hence,  in  the  WT/Tg  and  WT/null  heterozygotes  explored 
in  Chesebro  et  al.  (3),  we  anticipate  in  each  case  the  mem¬ 
brane  bound  PrP^j.  concentration  to  be  about  half  that  of  the 
homozygous  WT  mice,  whereas  the  intercellular  PrPjg 
concentration  should  be  about  half  that  of  the  homozygous 
Tg  mice.  Because,  as  shown  with  infectious  prions  bound 
to  electrodes  (7),  templating  and  conversion  can  be  driven 
by  scrapies  material  not  bound  to  the  membrane  surface,  we 
expect  the  incubation  time  of  the  WT/Tg  heterozygotes  to  be 
significantly  accelerated  relative  to  the  WT/null  heterozy¬ 
gotes  as  is  observed. 

We  note  that  membrane  associated  exponential  growth 
might  be  due  to:  i),  as  yet  undiscovered  membrane  specific 
enzymes  splitting  aggregates,  in  analogy  to  the  role  of 
Hp  104a  for  yeast  prions  ( 1 );  ii),  mechanical  breakage  of  aggre¬ 
gates  due  to  membrane  curvature  or  membrane  undulations 
(8);  iii),  oligomeric  autocatalysis  arising  from  interneuronal 


Biophysical  Journal:  Biophysical  Letters 


L78 


33 


Biophysical  Journal:  Biophysical  Letters 


templating  by  infectious  oligomeric  seeds  bound  to  one  or  the 
other  membrane  (9). 
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Abstract 

Numerous  human  diseases  are  associated  with  conformational  change  and 
aggregation  of  proteins,  including  Alzheimer’s,  Parkinson’s,  prion  diseases  (such  as  mad 
cow  disease),  familial  amyotrophic  lateral  sclerosis  (ALS,  or  Lou  Gehrig’s  disease), 
Huntington’s,  and  type  II  (mature  onset)  diabetes.  In  many  cases,  it  has  been 
demonstrated  that  conformational  change  and  aggregation  can  occur  outside  living  cells 
and  complex  biochemical  networks.  Hence,  approaches  from  materials  and  physical 
science  have  enhanced  our  understanding  of  the  role  of  protein  aggregation  in  these 
diseases  at  the  molecular  and  nanoscale  levels.  In  this  article,  we  will  review  what  is 
known  about  these  protein  structures  from  the  perspective  of  materials  science,  focusing 
on  the  details  of  emergent  oligomeric  and  nanotube-like  structures,  their  interactions  with 
model  lipid  bilayers,  how  the  structures  relate  to  observed  biological  phenomena,  and 
how  protein  aggregation  and  amyloid  formation  can  be  employed  for  the  good  in  biology 
and  materials  science. 

Keywords:  amyloid  diseases,  complex  adaptive  matter,  emergent  behavior,  prion 
diseases,  protein  aggregation,  protofibrils,  nanotubes,  nanowires. 


Amyloid  Diseases:  A  Truly 
Emergent  Phenomenon 

The  primary  goal  of  the  study  of  com¬ 
plex  adaptive  matter  is  to  identify  the  key 
organizing  principles  that  govern  mate¬ 
rials  phenomena  at  given  length  and  time 
scales.  In  the  case  of  the  folding  of  indi¬ 
vidual  globular  proteins  into  their  func¬ 
tional  structures,  we  know  a  great  deal. 
For  example,  the  local  rules  governing  a- 
helical  secondary  structure  have  been  un¬ 
derstood  in  detail  for  some  40  years;1,2 
indeed,  the  a-helix  was  predicted  by  Paul¬ 
ing  and  collaborators  prior  to  discovery.3 
For  the  global  tertiary  structure  of  the  pro¬ 
tein,  extensive  experiments  and  simula¬ 
tion  studies  show  that  proteins  engineered 
by  evolution  to  experience  minimal  frus¬ 
tration  in  the  interactions  between  closely 
contacted  amino  acids  (or  residues)  de¬ 
velop  funneled  energy  landscapes  and  rel¬ 
atively  rapid  (and  multiple)  pathways  to 
folding  after  synthesis.4  As  discussed  in 
the  article  by  Ramirez  in  this  issue,  "frus¬ 
tration"  refers  to  the  effect  of  competing 
interactions  that  make  it  impossible  to  fa¬ 
vorably  lower  the  interaction  energy. 


However,  we  have  comparatively  little 
understanding  of  the  organizing  princi¬ 
ples  governing  structure  formation  for 
proteins  interacting  with  other  proteins  or 
membranes.  This  lack  of  knowledge  is 
problematic,  because  proteins  left  alone 
tend  to  spontaneously  aggregate,  often  by 
formation  of  (3-sheet  structures,  which  are 
not  governed  by  local  formation  rules  like 
a-helices  (the  distances  between  hydro¬ 
gen  bonding  residues  along  the  backbone 
can  be  great).  |3-sheets  are  formed  from 
approximately  linear  stretches  of  peptide 
that  hydrogen-bond  from  line  to  line. 
These  structures  are  especially  prone  to 
protein  aggregation  due  to  favorable 
edge-to-edge  hydrogen  bonding  between 
sheets.5  This  aggregation  tendency  is  an 
obstacle  in  high-throughput  proteomics, 
where  one  is  interested  in  measuring  the 
properties  of  individual  proteins.6 

Purposeful  biological  aggregation  of 
monomeric  proteins,  as  in  the  assembly  of 
actin  filaments  or  microtubules,7  is  usually 
highly  regulated  and  energetically  con¬ 


trolled  (we  will  discuss  tightly  controlled 
and  biologically  useful  (3-sheet  self- 
assembly  later  in  the  article).  Living 
organisms  have  evolved  an  effective 
quality-control  system  to  prevent  protein 
misfolding  and  aggregation,  where  chap¬ 
erone  proteins  provide  "safe  houses"  for 
folding  proteins,8  and  the  ubiquitin/pro- 
teosome  system  ensures  rapid  degrada¬ 
tion  or  disposal  of  misfolded  proteins. 

Table  I  conveys  the  tragic  side  of  uncon¬ 
trolled  (3-sheet  self-assembly:  it  summa¬ 
rizes  the  key  aspects  of  seven  (out  of 
dozens)  prominent  human  amyloid  dis¬ 
eases.  Amyloid  means  "starch-like" — the 
aggregates  stain  like  starch.  These  dis¬ 
eases  typically  arise  in  old-  or  middle- 
aged  populations,  and  frequently  arise 
spontaneously  or  sporadically  rather  than 
from  genetic  predisposition.  Indeed,  for 
spontaneous  Alzheimer's,  Parkinson's, 
prion,  and  immunoglobulin  light  chain 
diseases,  incidence  varies  little  between 
countries.  Moreover,  prion  diseases,  the 
lone  infectious  type  of  amyloid  disease, 
show  highly  reproducible  dose-versus- 
incubation-time  distributions  for  inter- 
cerebrally  inoculated  animals.9  Remarkably, 
infectious  protein-only  prion  aggregates 
have  been  grown  in  vitro, 10  proving  that 
the  protein  interaction  properties  alone, 
without  additional  biochemical  guidance, 
account  for  the  disease. 

These  observations  suggest  that  these 
diseases  can  be  studied  from  the  perspec¬ 
tive  of  materials  growth,  without  exten¬ 
sive  biological  modulation.  Indeed,  it 
appears  that  the  growth  of  amyloid  struc¬ 
ture  has  much  in  common  with  the  ori¬ 
ented  aggregation  of  inorganic,  nearly 
monodispersed  nanoparticles.11  We  seek 
here  to  portray  the  growing  scientific 
movement  toward  the  use  of  concepts  and 
tools  from  materials  science  in  the  study  of 
amyloidogenic  proteins  to  elucidate  the 
mechanisms  of  disease  and  design  new 
materials. 

Amyloid  Structures:  Plaques, 
Protein  Nanotubes,  and  Oligomers 
Plaques 

The  extracellular  and/or  intracellular 
accumulation  of  amyloid  fibrils  in  the 
form  of  plaques  or  inclusions  (Lewy  bod¬ 
ies)  in  the  brain  is  a  defining  hallmark  for 
several  neurodegenerative  diseases  (Fig¬ 
ure  1).  For  example,  Alzheimer's  disease 
patients  have  large  quantities  of  postmortem 
brain  plaques,  predominantly  composed 
of  40-42-amino-acid-long  Ap  peptides 
that  are  cleaved  by  protease  proteins  from 
the  Alzheimer's  precursor  protein  (BAPP).13 
These  micrometer-scale  plaques  are  com¬ 
posed  of  multi-polymeric  strands  of  the 
peptide,  called  fibrils  (Figure  le),  which 
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Table  I:  Aspects  of  Some  Human  Amyloid  Diseases  with  Associated  Protein  Aggregation. 


Disease 

Alzheimer’s51 

Protein/Peptide 

p-42  (from  BAPP)  and  x 

Function 

? 

Heritable 

Component 

~25% 

Incidence 

—50%  of  post-85-yr-old 
population 

Onset  Age* 
(Years) 

>65 

Parkinson’s52 

a-synuclein 

? 

5-10% 

-1%  of  post-50-yr-old 
population 

55-60  yrs 

Huntington’s23 

huntingtin 

? 

100% 

1  in  20,000  (Caucasian) 

35-40 

Familial  ALS 
(Lou  Gehrig’s 
disease)53 

Superoxide  dismutase 
(20%  of  cases) 

Lowers  oxidative 

stress 

100%  (5-10%  of 
all  ALS) 

2  in  106 

46 

Type  II  diabetes54 

IAPP 

? 

High  (obesity 
trigger) 

14  x  106per  year  (U.S.) 

>40 

Immunoglobulin 
light  chain55 

IG  light  chain 

Immune  response 

Small/unknown 

1  in  1 0s 

64 

Prion  diseases35 

Prpc 

Lowers  oxidative 

stress? 

10-15% 

-1  in  106 

63 

*lncidence/onset  age  is  for  non-heritable  sporadic  disease,  unless  otherwise  noted. 


have  high  quantifies  of  P-sheet  structure 
as  revealed  by  crystallography,  circular 
dichroism  (which  detects  the  different 
light  polarization  rotation  tendencies  of  a- 
helices  and  P-sheets),  and  more  recently, 
solid-state  magnetic  resonance  studies. 
These  P-strands  are  aligned  perpendicular 
to  the  fibrillar  axis,  in  a  so-called  "cross-P" 
structure,  shown  in  Figure  If.  Postmortem 
plaques  and  inclusions  from  a  variety  of 


diseases  are  shown  in  Figure  l.13  We  note 
that  the  amyloid  fibrils  are  protein  nano¬ 
tubes,  hollow  in  the  middle,  with  diame¬ 
ters  of  the  order  of  10-20  nm. 

Despite  a  dominant  research  focus  on 
both  structure  measurements  and  com¬ 
puter  modeling  on  plaques,  an  emerging 
perspective  is  that  the  plaques  may  repre¬ 
sent  disease  end  points  having  little  to  do 
with  toxicity.14  This  view  is  supported  by 


observations  such  as  (1)  the  abundant  AP 
plaques  observed  in  the  brains  of  individ¬ 
uals  displaying  no  symptoms  of  Alzheimer's 
disease  and  (2)  the  non-uniformity  of 
plaque  observation  in  prion  diseases:  for  ex¬ 
ample,  victims  of  kuru  (a  disease  among 
the  Fore  people  of  New  Guinea,  arising 
from  ritual  cannibalism  of  deceased  tribe 
members)  exhibit  them,  while  victims  of 
spontaneous  Creutzfeldt-Jakob  disease 
usually  do  not.15  In  this  view,  toxicity 
is  engendered  by  small  p-sheet  aggre¬ 
gates,  possibly  on  the  pathway  to  amyloid 
fibrils. 

P-, Helices 

Accordingly,  attention  has  turned  to 
small  P-sheet  motifs  with  multiple  assem¬ 
bly  outcomes,  especially  the  left-handed 
p-helix  (LHBH)  structures  shown  in  Fig¬ 
ure  2, 16  recently  proposed  as  the  P-sheet 
unit  for  infectious  mammalian  prion 
trimers  on  the  basis  of  cryogenic  electron 
microscopic  data.17  The  LHBH  P-sheet 
structure  has  also  been  proposed  for  yeast 
prion-like  proteins,18  Alzheimer's  disease,19 
and  Huntington's  disease.16'20  This  LHBH 
motif  was  first  observed  in  several  bacter¬ 
ial  enzymes  and  the  "antifreeze"  protein 
of  the  spruce  budworm;  to  date,  there  are 
11  structures  in  the  protein  data  bank 
(http://www.rcsb.org/pdb/)  confirmed 
to  have  LHBHs.  LHBHs  are  usually  pre¬ 
sumed  (or,  in  two  cases,  observed)21'22  to  be 
in  protein  trimers.  The  LHBH  motif  has  a 
fundamental  repeat  unit  of  triangular 
cross  section,  consisting  of  18  amino  acids 
with  two  per  bend  region  and  four  per 


Figure  1.  Plaques  and  fibrils,  (a)-(d)  Postmortem  tissue  plaques  from  human  amyloid 
disease  patients  (from  Reference  13):  (a)  Ap42  (Alzheimer's)  plaque,  (b)  a-synuclein 
(Parkinson's)  plaque,  (c)  huntingtin  plaque  (Huntington's  disease),  and  (d)  PrPSc  (kuru) 
plaque,  (e)  Transmission  electron  microscope  images  of  amyloid  fibrils  (H.  Lashuel, 
unpublished  data),  (f)  Model  of  a  hollow-core  SH3  domain  fibril:  upper  panel  shows  a 
density  map  cross  section,  lower  panel  shows  the  cross-p  structure.56 
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Figure  2.  (a)  Peptide  backbone  of  measured  left-handed  p-helix  (LHBH)  structure  of  1 FWY 
protein,  (b)  Proposed  LFiBFi  structure  of  infectious  prion  protein  (after  Reference  16). 
Parallel,  hydrogen-bonded  p-sheets  are  shown  as  strips,  (c)  Superposed  structure  of  six 
LHBFI  backbones,  illustrating  the  uniform  cross  section  of  the  helical  structure  (red,  1  FWY 
protein;  blue,  proposed  infectious  prion  protein;  yellow,  1 HMO  protein;  white,  1J2Z  protein; 
green,  1L0S  protein;  fuchsia,  1 T3D  protein),  (d)  LHBFI  “nanotube”  formed  from  a  stack  of 
three  1FWY  helices.  All  images  generated  with  the  RASTOP  molecular  viewing  program 
(see  www.geneinfinity.org/rastop/);  all  files  for  the  proteins  1FWY,  1HMO,  1J2Z,  1  LOS,  and 
1 T3D  input  to  RASTOP  are  from  the  Protein  Data  Bank  (http://www.rcsb.org),  except  for 
the  PrPSc  model,  which  is  courtesy  of  C.  Govaerts. 


(3-strand,  alternating  between  hydropho¬ 
bic  and  hydrophilic  in  the  strand.  The 
triangular  cross  section  is  very  uniform 
across  observed  LHBHs  (Figure  2c); 
each  edge  is  ~1.9  nm  in  length.  In  Fig¬ 
ure  2d,  we  show  how  multiple  copies 
of  the  LHBH  of  an  enzyme  from  E.  coli 
(labeled  1FWY  on  the  protein  data  base 
at  www.rcsb.org/pdb)  can  be  assembled 
into  a  kind  of  (3-nanotube.  This  raises 
the  natural  and  intriguing  question: 
are  amyloid  fibrils  composed  of  such 
"nanofilaments"?19 

The  LHBH  has  a  tantalizing  connection 
to  Huntington's  and  other  polyglutamine 
repeat  diseases,  where  pathology  derives 
from  an  inherited  excess  number  of  re¬ 
peats  of  the  amino  acid  glutamine  on  one 
end  (the  N-terminus)  of  the  huntingtin 
protein.23  A  polyglutamine  repeat  number 
p  of  less  than  ~24  is  normal;  for  p  >  36, 
disease  is  certain.  Note  that  36  is  the  num¬ 
ber  of  amino  acids  or  residues  in  two  turns 
of  a  LHBH  that  has  fully  saturated  inter¬ 
nal  hydrogen  bonds.19 

The  LHBH  presents  a  challenge  to 
theory  and  simulation.  At  present,  this 
motif  has  not  emerged  from  any  molecu¬ 
lar  dynamics  simulations  or  semi-analytic 
(Hamiltonian  or  cellular  automata)  ap¬ 


proaches.  The  complicating  factor  is  the 
long-range  coupling  along  the  backbone 
(bonded  amino  acids  are  separated  by  18 
residues).  On  the  other  hand,  the  remark¬ 
able  conservation  of  the  shape  and  helical 
cross  section  suggest  that  this  is  a  motif 
ripe  for  study. 

Nanoscale  Oligomers 

(3-sheet  oligomers  obtained  from  in  vitro 
growth  of  aggregates  have  been  studied 
extensively.13  One  common  form  seems  to 
be  a  spherical  micelle-like  aggregate  of  (3- 
converted  proteins,  which  have  been  im¬ 
plicated  as  precursors  for  the  chain-like 
(protofibril)  and  annular  oligomers  that 
are  also  seen  frequently  during  in  vitro 
amyloid  formation  by  most  amyloido- 
genic  proteins.24  Figure  3  shows  a  com¬ 
pendium  of  oligomers  and  protofibrils 
from  amyloid  disease  proteins.25,26  The  an¬ 
nular  oligomers  underlie  the  proposed 
toxicity  mechanism  discussed  later  in  this 
article.  We  note  that  all  atomic  force  micro¬ 
scope  (AFM)  measurements  of  spherical, 
chain-like,  and  annular  oligomers  share 
the  same  height  and  diameter,  suggesting 
that  the  spherical  aggregates  are  the  pre¬ 
cursors  to  the  larger  chain  and  annular 
pore-like  structures. 


The  diversity  of  pre-fibrillar  oligomer 
structures,  the  formation  of  which  are 
highly  dependent  upon  protein  sequence 
and  environmental  conditions  (e.g.,  pH, 
salt  concentration,  and  levels  of  molecular 
crowding),  may  explain  the  mystery  of 
prion  disease  strains.  Strains,  for  a  given 
mammal,  have  unique  incubation-time 
versus  dose  distributions,  tissue  lesion  pro¬ 
files,  and  distributions  of  post-translationally 
attached  sugars:  the  prion  protein  can 
have  0,  1,  or  2  sugars  attached.27,28  More¬ 
over,  strains  breed  true  upon  multiple 
passage  in  animals.  This  means  that  upon 
passage  from  a  diseased  animal  to  a 
healthy  animal,  the  same  properties  (incu¬ 
bation  time,  lesion  profile,  sugar  binding) 
are  preserved.  There  is  considerable  evi¬ 
dence  that  strain  is  encoded  in  prion  con¬ 
formation  but  no  detailed  understanding 
of  the  underlying  mechanisms.27,28  Prion 
oligomers17  might  have  a  spectrum  of  dif¬ 
ferent  shapes,  and  oriented  aggregation  of 
such  protein  "nanoparticles"  might  "breed" 
the  conformation  true  (nonmatching  shapes 
will  be  energetically  unfavorable). 

Domain  Swapping 

A  separate  way  to  generate  fibrils  and 
oligomers  is  by  "domain  swapping,"29,30 
in  which,  say,  two  identical  copies  (A,B)  of 
a  protein  exchange  a  domain  (a  well-defined 
protein  region  attached  to  a  flexible  sec¬ 
tion  of  peptide).  The  domain  of  monomer 
A  binds  to  the  corresponding  region  of  B 
and  vice  versa.  The  swapping  is  not  lim¬ 
ited  to  dimers;  it  can  lead  to  filamentary 
structures  in  which  the  zth  protein  swaps 
with  the  (i  +  l)th  protein,  for  example, 
closed-chain  oligomers  (where  filament 
ends  are  brought  together  to  domain- 
swap)  or  two-dimensional  structures. 

Domain  swapping  has  been  proposed 
to  play  a  role  in  prion  strains.31  Domain 
swapping  in  protein-protein  interactions 
generically  and  amyloid  diseases  in  par¬ 
ticular  seems  certain  to  emerge  as  a  critical 
theme  in  the  coming  years.  At  the  concep¬ 
tual  level,  the  study  of  coarse-grained  pro¬ 
tein  models  with  molecular  dynamics  has 
shown  that  the  formation  of  domain- 
swapped  dimers  will  proceed  down  a 
funneled  landscape  if  the  dimer  enjoys 
minimal  frustration. 

Aggregation  Pathways  and 
Kinetics 

AFM  and  kinetic  modeling,  well  known 
to  materials  scientists,  have  been  instru¬ 
mental  in  advancing  our  understanding 
of  the  structural  properties  of  the  protein 
aggregates  linked  to  disease  and  their 
growth  kinetics.  Figure  4  schematically 
shows  accepted  models  on  amyloid  con¬ 
version/aggregation  kinetics.  We  note 
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Figure  3.  Amyloid  oligomers,  (a)  Ap  peptide  implicated  in  Alzheimer's  disease,  observed  by 
electron  microscopy  (scale,  50  nm/box)  (from  Reference  24).  (b)  Nine  different  oligomers  of 
mutant  a-synuclein,  implicated  in  Parkinson's  disease,  observed  by  electron  microscopy 
(scale,  40.5  nm/box)  (from  Reference  25).  (c)  Six  oligomers  of  mutant  SOD1,  implicated  in 
familial  ALS  (from  Reference  56).57  (d)  Nine  mutant  a-synuclein  oligomers  (scale,  40.5  nm/box) 
(from  Reference  25).  (e)  Atomic  force  microscope  (AFM)  images  ofAp  oligomers  inserted 
in  a  supported  lipid  bilayer  (from  Reference  26).  Clearly  resolved  annular  oligomers  have 
outer  diameters  of —16  nm.  (1)  AFM  image  of  a-synuclein  oligomer  on  supported  lipid 
bilayer  (from  Reference  25).  (g)  Schematic  illustration  of  oligomerization  pathways.  Flere, 
the  subscript  n  in  the  second  stage  represents  the  number  of  monomers  in  an  oligomer. 


that  misfolds  in  the  diseased  form.37 
Given  a  mean  incubation  time  for 
kuru  of  12  years,38  this  suggests  that  en¬ 
demic  sporadic  prion  disease  requires 
~  1000-year  life  spans!  Meanwhile,  the  1  in 
106  sporadic  disease  background  inci¬ 
dence  may  reflect  the  low-amplitude, 
pre-peak  tail  in  the  incubation  time 
distribution.38 

Amyloid-Membrane  Interaction 
and  Toxicity 

Many  amyloidogenic  proteins  associate 
with  lipid  membranes.  AFM  studies  on 
supported  bilayers  and  molecular  model¬ 
ing  have  helped  support  a  potential  unify¬ 
ing  hypothesis  for  amyloid  disease 
toxicity:  that  small  oligomers  pierce  cell 
membranes,  triggering  cell  death  through 
superfluous  ion  pores.  It  was  proposed  in 
detail  that  a  pair  of  P-coupled  annular 
tetramers  of  Ap  peptides  can  insert  into 
the  neuronal  membrane  leaflets  and  create 
an  ion  pore.39  In  vitro- grown  oligomers 
qualitatively  consistent  with  this  hypothe¬ 
sis  have  been  found  (Figure  3)  and  studied 
with  AFM  on  supported  bilayers  (Fig¬ 
ure  3f);  these  are  evidently  composed  of 
4  nm  spherical  oligomers,  much  larger 
than  the  peptide  tetramers  of  Refer¬ 
ence  39.  These  oligomers  permit  excess 
calcium  flux,  which  is  toxic  to  cultured 
neurons.26  Intracerebral  inoculation  of  rats 
and  mice  with  a  solution  presumed  rich  in 
oligomeric  Ap  particles  leads  to  a  re¬ 
versible  short-term  memory  deficit.40  A 
recent  theoretical  study  found  strong  cor¬ 
relations  between  the  membrane  insertion 
configuration  and  the  pore  model  for  4 
of  5  mutations  leading  to  early-onset 
Alzheimer's  disease.41  While  the  pore 
model  is  not  universally  accepted,  and 
may  not  apply  to  all  the  diseases  (the 
prion  trimer  model,  for  example,  will  not 
allow  ion  passage),  it  remains  a  vibrant 
area  of  research. 


that  P-sheet  conversion  is  typically  con¬ 
comitant  with  aggregation  (i.e.,  conforma¬ 
tional  changes  and  protein  assembly  are 
linked),32-33  with  the  possible  exception  of 
polyglutamine  proteins,34  and  requires 
templating  either  by  spontaneously  formed 
(and  rare)  nuclei  or  by  external  seeding  of 
aggregates.33  Also,  fission  of  mammalian 
prions  is  necessary  for  disease  propagation 
and  exponentially  growing  aggregates  in 
mammals;  prion-like  proteins  in  yeast 
have  fission  effected  by  another  protein35 
(although  the  fission  can  occur  sponta¬ 
neously  in  vitro).36  Why  mammalian  prions 
fission  and  other  amyloidogenic  proteins 
do  not  remains  a  mystery. 

In  vitro  and  theoretical  studies  of  con¬ 
version/aggregation  suggest  that  the  late- 


age  onset  in  amyloid  diseases  derives  from 
slow  underlying  molecular  processes.  For 
example,  in  vitro  kinetics  experiments  for 
polyglutamine  peptides  extrapolated  to  in 
vivo  concentrations  of  huntingtin  protein 
suggest  that  within  the  "sampling  win¬ 
dow"  of  a  human  lifespan,  toxic  aggregate 
concentrations  should  arise  only  for  gluta¬ 
mine  number  p  >  36,  consistent  with  clin¬ 
ical  observations.34 

For  prion  diseases,  theoretical  modeling 
of  two-dimensional  aggregation  and  fis¬ 
sion  (prions  live  mostly  on  neuronal 
membranes)  yielded  a  sporadic  incuba¬ 
tion  time  distribution  that  peaked  at  ~  100 
times  that  obtained  from  dilute  seeding 
for  physiological  concentrations  of  the 
normally  expressed  protein  called  PrPc 


Amyloids  for  Good  in  Biology  and 
Materials  Science 
Biologically  Useful  Amyloids 

Emerging  evidence  suggests  that  amy¬ 
loid  structures  can  provide  useful  biologi¬ 
cal  functions.  Some  examples  are 

1.  Heritable  amyloid  structure  in  yeast.  As 
alluded  to  previously,  prion-like  proteins 
in  yeast  form  aggregates  that  fission 
upon  cell  division,  can  actively  confer 
phenotype,  and  may  provide  some  stress 
protection.42 

2.  Spider  silk.  Spiders  produce  insoluble  fil¬ 
aments  of  fibroin  protein  that  possess  30x 
greater  extensibility  and  toughness  than 
steel.  Recent  circular  dichroism  studies 
show  that  significant  amyloid-like  cross-P 
structure  develops  in  a  region  of  reduced 
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Figure  4.  Schematic  illustration  of  amyloid  conversion  and  aggregation.  Native  monomers  of 
the  relevant  peptides  or  proteins  are  shown  as  circles  and  high-p-content  converted 
proteins  as  squares.  Peptides  can  either  spontaneously  convert  or  be  seeded.  Subsequent 
aggregation  and  conversion  can  generate  either  proto-fibrils  via  the  elongation  step,  which 
then  form  amyloid  fibrils  and  plaques,  or  pass  to  “off-pathway”  oligomers.  In  the  case  of 
prions,  the  oligomers  and/or  fibrils  can  fission,  which  then  re-seeds  the  initial 
conversion/aggregation  process.  Prions  are  the  only  known  amyloid  proteins  to 
spontaneously  fission  in  the  converted  form  either  in  vitro  or  in  vivo.  It  is  assumed  (but  not 
shown  explicitly  here)  that  both  oligomers  and  native  monomers  experience  some  level  of 
removal  by  cellular  systems  in  vivo. 


pH  downstream  from  the  initial  extrusion 
site.43 

3.  Chorion  in  fish  and  insect  egg  shells.  Trun¬ 
cated  peptides  from  the  central  regions  of 
two  chorion  proteins  self-assembled  into 
spherulites  possessing  (3-structure,  which 
then  converted  to  fibrils  upon  maturation, 
suggesting  amyloid  character  in  actual  egg 
shells.44 


4.  Amyloid-like  structure  in  synapses  of 
Aplysia  (slugs).  The  N-terminus  region  of 
the  Aplysia  synapse  protein  CPEB  is 
glutamine-rich,  like  huntingtin  protein. 
Engineered  expression  of  CPEB  in  yeast 
yields  prion-like  aggregates  similar  to  the 
native  ones  discussed  earlier;  hence,  prion¬ 
like  states  of  CPEB  might  effect  long-term 
strengthening  of  synaptic  contacts.45 


Amyloids  in  Materials  Science 

The  regular  diameter  and  periodicity  of 
the  amyloid  fibrils  make  them  good  mate¬ 
rials  templates.  One  group  employed  yeast 
prion-like  proteins  to  template  ~100-nm- 
wide  gold  nanowire  growth;  the  gold- 
coated  protein  filaments  after  initial  decor¬ 
ation  by  small  gold  nanoparticles  were 
found  to  bind  to  genetically  engineered 
cysteine  residues.46  Silver  nano  wires  of 
20  nm  width  were  grown  inside  filamen¬ 
tary  cross-(3  tubes  grown  from  a  dipheny¬ 
lalanine  peptide.47  Engineering  of  con- 
trollably  switched  (3-sheet  materials  could 
prove  valuable  for  tissue  growth  scaffold¬ 
ing,  as  one  example.48  Hybrid  molecules 
including  eight  amino  acids — four  per 
strand,  but  with  non-amino  bends — have 
been  developed  that  controllably  self- 
assemble  into  different  (3-sheet  structures 
(as  shown  in  Figure  5)  depending  upon 
the  pH,  analogous  to  spider  silk.49  Clearly, 
the  future  of  engineered  amyloid  struc¬ 
tures  in  materials  science  looks  bright.50 

Conclusion 

In  this  article,  we  have  developed 
themes  that  bring  together  the  fields  of 
amyloid  diseases  and  materials  science. 
Protein  misfolding  and  aggregation  phe¬ 
nomena  are  intimately  linked  to  many  se¬ 
rious  public  health  issues.  However,  many 
aspects  of  the  phenomena  have  close 
analogies  in  synthesized  materials,  and 
their  full  molecular  understanding  re¬ 
quires  experimental  and  modeling  tools 
more  familiar  in  the  physical  and  mate¬ 
rials  sciences.  In  addition  to  the  possible 
medical  breakthroughs  that  such  multi¬ 
disciplinary  studies  can  lead  to,  there  is 
a  growing  possibility  that  understand¬ 
ing  the  mechanisms  of  amyloid  formation 
can  have  wide  impact  in  fields  ranging 
from  basic  neurobiology  to  materials 
science. 


Figure  5.  (a)-(c)  Images  from  electron  microscopy  showing  pH-dependent  fibril  growth  of  designed  peptidomimetic  molecules  (from 
Reference  49).  Scale  bars  are  100  nm. 
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Sort  Results  by: 

1.  California  researcher  finds  that  mad  cow  infection  spreads  rapidly 

Scottsbluff  Star-Herald  -  Jun  1 6  8:06  AM  ^ 

DAVIS,  Calif.  (AP)  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow 
so  frighteningly  fast. 

2.  Modeling  to  Understand  Mad  Cow  Disease 

SAP  Info -Jun  16  2:15  AM  ^ 

Two  physics  professors  at  the  University  of  California  Davis,  who  have  been  researching  mad  cow  disea 
use  computerized  mathematical  models  to  explain  how  the  disease  acts  so  fast. 

3.  Researchers  explore  spread  of  mad  cow  disease 

CTV.ca  -  Jun  15  4:49  PM  ^ 

DAVIS,  Calif.  —  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  dis 
frighteningly  fast  using  an  unconventional  tool  mathematical  modeling. 

4.  UC  Researchers  Explore  Spread  of  Mad  Cow 

TechNewsWorld.com  -  Jun  15  8:27  AM  ^ 

Scientists  believe  the  human  body  usually  rejects  misfolded  proteins  without  harm;  it's  when  a  microscopi 
flawed  proteins  is  consumed  from  another  animal  that  the  disease  can  take  hold  in  a  person  or  animal  wh 
would  never  develop  the  disease. 

5.  The  Capital  City's  Newspaper  Online 
Independent  Record  -  Jun  14  11 :1 2  PM  ^ 

AP  photo  -  Rajiv  Singh,  a  physics  professor  at  the  University  of  California  Davis,  holds  up  a  slide  showint 
of  a  human  protein  cell  at  his  campus  office  in  Davis,  Calif.,  on  Tuesday. 

6.  Mad  cow  computer  modeling  finds  infection  spreads  rapidly 

Billings  Gazette  -  Jun  15  7:31  AM  ^ 

DAVIS,  Calif.  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  dise 
frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

7.  Mad  cow  computer  modeling  finds  infection  spreads  rapidly 

Billings  Gazette  -  Jun  15  7:32  AM  ^ 

DAVIS,  Calif.  --  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  dise 
frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

8.  UC  Researchers  Explore  Spread  of  Mad  Cow 

First  Coast  News  -  Jun  15  10:04  AM  ^ 

DAVIS,  CA  (AP)  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  < 
frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

9.  Mad  cow  computer  modeling  finds  infection  spreads  rapidly 
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SanLuisObispo.com  -  Jun  15  3:35  AM 

DAVIS,  Calif.  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  dise 
frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

10.  UC  researchers  explore  spread  of  mad  cow  ^ 

Kansas  City  Star  -  Jun  1 5  5:20  AM 

DAVIS,  Calif.  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  dise 
frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

11.  UC  Researchers  Explore  Spread  of  Mad  Cow ^ 

RedNova  -  Jun  15  7:26  AM 

DAVIS,  Calif.  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  dise 
frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

12.  UC  researchers  explore  spread  of  mad  cow  ^ 

The  Times  Leader  -  Jun  15  5:18  AM 

DAVIS,  Calif.  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  dise 
frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

13.  UC  researchers  explore  spread  of  mad  cow  ^ 

Macon  Telegraph  -  Jun  15  5:28  AM 

DAVIS,  Calif.  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  dise 
frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

14.  Mad  cow  computer  modeling  finds  infection  spreads  rapidly  ^ 

Monterey  County  Herald  -  Jun  15  3:29  AM 

DAVIS,  Calif.  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  dise 
frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

15.  UC  researchers  explore  spread  of  mad  cow  ^ 

TimesLeader.com  -  Jun  15  5:08  AM 

DAVIS,  Calif.  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  dise 
frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

16.  Mad  cow  computer  modeling  finds  infection  spreads  rapidly  ^ 

The  Times  Leader  -  Technology  -  Jun  15  5:48  AM 

DAVIS,  Calif.  (AP)  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow 
so  frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

17.  UC  Researchers  Explore  Spread  of  Mad  Cow  ^ 

WJLA-TV  Washington  D.C.  -  Jun  15  6:08  AM 

DAVIS,  Calif.  (AP)  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow 
(  news  )  acts  so  frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 

18.  UC  researchers  explore  spread  of  mad  cow  ^ 

Kansas  City  Star  -  Jun  15  5:27  AM 

DAVIS,  Calif.  -  Researchers  at  the  University  of  California,  Davis  are  trying  to  explain  how  mad  cow  dise 
frighteningly  fast  using  an  unconventional  tool  -  mathematical  modeling. 
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